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Abstract 

In this paper we prove the non-Unear asymptotic stability of the five- 
dimensional Schwarzschild metric under biaxial vacuum perturbations. This 
is the statement that the evolution of {SU (2) x U (l))-symmetric vacuum per- 
turbations of initial data for the five-dimensional Schwarzschild metric finally 
converges in a suitable sense to a member of the Schwarzschild family. It con- 
stitutes the first result proving the existence of non-stationary vacuum black 
holes arising from asymptotically flat initial data dynamically approaching a 
stationary solution. In fact, we show quantitative rates of approach. The proof 
relies on vectorfield multiplier estimates, which are used in conjunction with 
a bootstrap argument to establish polynomial decay rates for the radiation 
on the perturbed spacetime. Despite being applied here in a five-dimensional 
context, the techniques are quite robust and may admit applications to various 
four-dimensional stability problems. 
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1 Introduction 

The existence of black holes features among the most fundamental predictions of 
general relativity. In the appropriate mathematical language of the theory, these 
objects correspond to solutions of the Einstein equations 

1 

Rfiy - -^Rgfii' = SttT^j^ (1) 

possessing a regular event horizon and a complete null- infinity. General relativity 
admits an initial-value formulation suggesting that the appropriate setup to study 
black holes is in evolution from initial data. In this context, the main objective is 
to determine whether the maximal development associated to given data admits a 
complete null-infinity and a regular event horizon. 

Some important special black hole solutions (hence their initial-data) are known 
in closed form. They are static or stationary, with the well-known Schwarzschild 
and Kerr family of solutions amongst them, which are believed to play crucial roles 
as "final states" in gravitational collapse. It is fundamental for our understanding 
of the theory to investigate the stability of these explicit solutions, that is to say 
the global structure of the evolution arising from initial data close (in an appropri- 
ate sense) to that of the known reference solution. Due to the complexity of this 
non-linear problem, most rigorous studies have been focussed on special symmetry 
classes. Specifically, a paramount problem of black hole physics, the full non-linear 
stability of the Kerr-solution, remains open to date. 

A model in which both the global spacetime structure associated to the evolu- 
tion of general initial data and the stability of certain solutions in particular have 
been mathematically understood previously is that of the self-gravitating scalar 
field under spherical symmetry. The assumption of spherical symmetry casts the 
Einstein equations as a l-fl dimensional system of PDEs, the inclusion of a mass- 
less scalar field being the simplest way to circumvent Birkhoff's theorem^ In the 
context of this model, Christodoulou [2] proved that generic initial data either dis- 
perse, i.e. asymptote to Minkowski space for late times, or collapse to regular black 
holes. His seminal work was extended by Dafermos and Rodnianski [6], who proved 
that the development of initial data collapsing to black holes in fact approaches a 
Schwarzschild-metric on the exterior of the black hole at a sufficiently fast polyno- 
mial rate. These decay rates [6] of the scalar field were first suggested on a heuristic 
level by Price [8], and are thought to be sharp. It is remarkable that [6] is a "large 
data" result. The initial data need not be assumed close to Schwarzschildean; all 
initial data containing a trapped surface are shown to approach a Schwarzschild 
metric. 

1.1 The model 

An alternative model allowing the study of gravitational collapse in vacuo was re- 
cently proposed by Bizon et al. [1]. To understand their idea we recall that, in 
view of the four-dimensional Birkhoff's theorem, gravitational collapse in vacuo 

^Birkhoff's theorem impUes that spherically symmetric vacuum solutions are either 
Minkowskian or Schwarzschildean. 
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(T^i, = in ([T])) cannot be studied under spherical symmetry. In axisymmetry 
on the other hand, the Einstein equations no longer reduce to a system of 1 + 1 
dimensional PDEs and the resulting problem does not seem tractable with cur- 
rent mathematical techniques. The way out of this dilemma suggested by p] is 
to study the Einstein vacuum equations under SU (2)-symmetry in five dimen- 
sions. This is motivated by the following observation: The analogue of spherical 
symmetry in four dimensions, i.e. an SO (3) action on an orbital two-sphere, is 
clearly an SO (4) = {SU (2)^ x SU (2)^) /Z^ action on a 3-sphere in five dimen- 
sions. However, via the latter isomorphism there exist subgroups of 5*0(4), for 
instance SU {2)^ and {SU {2)j^ x U (1)^) /Z^ which still act transitively on the 3- 
sphereH Consequently, even within the class of the smaller symmetry-groups (com- 
monly called triaxial- or biaxial- Bianchi IX depending on the subgroup to which 
one restricts) the Einstein equations reduce to a system of 1+1 dimensional PDEs. 
Moreover, Birkhoff's theorem is evaded by the introduction of one or two (in the 
triaxial case) dynamical degrees of freedom arising from the reduced symmetry. 

In the biaxial case this degree of freedom is manifest in a certain function _B, 
which geometrically speaking corresponds to the "squashing" of the three sphere. B 
is normalized such that it is zero for the Schwarzschild-Tangherlini metric. From the 
point of view of the analysis it can be understood as the analogue of the massless 
scalar field in four dimensions. The Einstein equations ([1]) imply the following 
non-linear wave equation for the squashing field B 



In [T] the model outlined was investigated numerically, suggesting that small 
initial data will disperse, whereas large data will collapse to black holes, approaching 
some Schwarzschild-Tangherlini black hole for large times. The mathematical study 
of the model was initiated shortly thereafter by M. Dafermos in collaboration with 
the present author. In [5], the following statemcnt[f| was proven: 

Theorem. Consider a triaxial- symmetric initial data set {Y,^g^K), which is close 
in an appropriate norr^ to an initial data set {"S, gs, Ks) evolving to the five di- 
mensional Schwarzschild-Tangherlini solution of mass M. Let the squashing fields 
Bi, B2 which are identically zero for the five- dimensional Schwarzschild metric, he 
of compact support on the initial hypersurface. Let Q he the Lorentzian quotient of 
the future Cauchy development of the data. Then Q contains a subset with Penrose 
diagram: 



It particular, the quotient of the maximal development of the set {Y,,g,K) admits 
a complete null-infinity with final Bondi mass M / close to M , and a regular event 
horizon 7i+ on which the Penrose inequality < 2A4 / holds. Here r is the area- 
radius function. 

The above theorem can be paraphrased as stating that perturbations of Schwarz- 
schild-Tangherlini initial data again collapse to regular black holes close to the 
original Schwarzschild black hole. This result was termed orhital stahility of the 

■^Thc subscripts L and R stand for the left and the right action respectively. 
^Actually, it follows from a stronger statement proven in [5]. 
*See [5] for the precise definition. 




(2) 
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fivc-dimcnsional Schwarzschild metric in [5] and generated the first vacuum black 
hole solutions arising from asymptotically flat initial data that arc not stationary!^ 
Crucial for the proof of the above theorem is the existence of good monotonicity 
properties for a function m{u,v), called the Hawking mass, defined in (|29p . It 
converges to the ADM mass defined at the asymptotically flat end. It is shown to 
satisfy 9„m < and d^m > on the domain of outer communications, leading to 
an a-priori bound for the total mass fluctuation on the spacetime in terms of the 
initial data. 

1.2 The main theorem 

Orbital stability provides of course certain control over the global structure of the 
solution. Nevertheless, it leaves the details of the late-time behaviour unclear. 
In particular, solutions could exhibit unexpected features at late times with the 
squashing field B oscillating in some complicated manner and the geometry thus 
never settling down. This problem is finally addressed in the present paper. By 
proving appropriate decay-rates we will show that the squashing field does decay 
for late times and hence that perturbations converge to another member of the 
Schwarzschild- Tangherlini family. 

1.2.1 The statement 

The main result is 

Theorem 1.1. Consider a biaxial- symmetric initial data set [Yi^g^K), which is 
close in the sense of the previous theorem to an initial data set {T,, gs, Ks) whose 
maximum development is the five dimensional Schwarzschild-Tangherlini solution 
of mass M . Let -k : Ai ^ Q denote the projection map of the maximal development 
of (E, g, K) to the two-dimensional Lorentzian quotient space Q and let S ^ tt (S). 
Fix a curve of constant area radius, r = rx, away from the horizon, intersecting S 
at P as depicted below. Assume furthermore that the initial data slice S coincides 




for r > rK with an integral curve of the globally defined vectorfield Vr on Q and 
that the data is Schwarzschildean outside a compact set, i.e. that the squashing field 
B is of compact support. 

Define regular coordinates {u,v) on the subset J+ (^S H {r > r/^}^ H i^^) 
of the Penrose diagram arising by the previous Theorem as follows. Let the point 
R, determined by the intersection of the curve = 4m with S, have coordinates 
u ~ V = \IM . Set r^y = i(l — /i), with fj, = along the null-ray PQ and 

r^u = along null-infinity. Ln these coordinates u ^ oo along null-infinity as 
i~^ is approached. The horizon TL^ is parametrized as {oo,v). Define t = '^^^ and 

„* v—u 

' 2 ■ 

^Solutions with a future complete, but not past complete, X+ have been constructed previously 
by Chrusciel [1], by solving a certain parabolic problem. 
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Then there exists a dimensionless constant S > 0, depending only on the geom- 
etry of S such that if the field B satisfies 



M-i 



r2 \B\ + r2 



7-2 



< 5 



on S f] {r > rK} and 
1 



M 



Sn{r>rK} 



{duBY + ^2 {d^BY + (u^ + v^) (-r,„) 



as well as 



B, 



< S 



r^\B\ + r2 

on the ray v ^ v (P) H {r < ta'} , then the squashing function B satisfies 



\B\+VM\B^,\ 
where v+ = max (1, f ), 

\b\ + Vm\bj 



CVM , 
< for r < rx 



B, 



CVM , 
< for r > rx 



(3) 



(4) 



(5) 



(6) 



(7) 



C Mi 

\B\ < 3 — for r>rK 



(8) 



on D ^ J+ (^"S n {r > n J (T^) for a dimensionless constant C (which de- 

pends on the choice of r^ ) computable from the initial data. 

Wc will refer to this result as the asymptotic stability of the Schwarzschild- 
Tangherlini solution. In particular, Theorem 11.11 produces the first dynam- 
ical vacuum solutions arising from asymptotically flat initial data and 
converging to stationary black holes for late times. 

1.2.2 Remarks 

Restricting 5* to coincide with a Vr integral curve for r > rif is justified by Cauchy 
stability and the fact that the global properties of the Penrose diagram are already 
known by the orbital stability result of the previous theorem. It has been assumed 
to avoid some clumsy notation in the proof. 

Cauchy stability also justifies stating the smallness assumptions Q, (21) and ([5]) 
on the slice 



Sr,, = (sc^{r> rx}) U {{v = v [P)} n {r < r^}) 



(9) 



instead of sE 

The advantage of doing it this way is that ^ and ([3 do not depend 
on the choice of double-null coordinates on the Penrose diagram|j Assumption (|4|) 
on the contrary depends on the choice of coordinates. However, since both the u 
and the v coordinate are easily shown to be finite in the region where B is supported 



®The smallness assumption ^ easily translates into an appropriate smallness assumption on 
S, depending on the geometry of S for r < rx, after extending the coordinate system to all of 

j+ {s)r\.J- (x+). 

^This will become useful later because the bootstrap argument applied in the proof requires 
the definition of different coordinate systems. 
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on S* n {r > vk} in the given coordinate system, assumption ([¥]) is automatically 
satisfied if we choose the (5 in small enough since B and is assumed to be of 
compact support initially. Hence it could be dropped by making 6 even smaller. We 
have nevertheless included ^ for conceptual reasons which will become apparent 
later in the proofH The condition ([4]) would also be required if one eventually drops 
the assumption of compact support, for then ([4]) imposes conditions on the decay 
of the fields near infinity. 

Factors of \/M have been inserted in all formulae to make constants dimension- 
less. 

1.3 Summary of the proof 

Before we embark upon an outline and a discussion of the proof, it is perhaps 
illuminating to compare and contrast the situation with the proof of Price's law 
[6] for a self-gravitating spherically symmetric scalar field in 3 + 1 dimensions. It 
turns out that the techniques developed in the latter paper to derive decay rates 
do not generalize to the system under consideration. The underlying reason can be 
traced back to two crucial estimates applied in [6]. The first of these, which allows 
one to extract decay directly from the horizon, relies heavily on the homogeneity of 
the non-linear wave equation satisfied by the field in the scalar field model. The 
second estimate is made possible by the existence of an almost Riemann invariant, 
a quantity admitting better decay properties than the scalar field </) itself, which can 
be exploited to derive uniform decay of the energy in the area radius r. This decay 
played an important role in conjunction with the pigeonhole principle completing 
the argument in [6]. 

In the five-dimensional case there is no almost Riemann invariant and hence no 
apparent analogue to obtain decay in r for the energy in the asymptotic region. 
Moreover, the wave equation ^ satisfied by the dynamical field B has an inhomo- 
geneous part, which in particular appears in the redshift estimate. These obstacles 
necessitate a very different approach to proving decay. The path we choose here 
is based on exploiting energy currents arising from vectorfield multipliers. This 
method was already central in the proof of the non-linear stability of Minkowski 
space [3] and has recently been applied at the linear level in the black hole context 
for the first time [7]. In the latter paper, decay rates for a scalar field satisfying 
the homogeneous linear wave equation on a four-dimensional Schwarzschild space- 
time are proven^ Key to establishing decay, at least away from the horizon, is the 
application of a so-called Morawetz vectorfield. A careful analysis reveals that the 
decay-rates can be generalized to the linear problem associated with the non-linear 
problem studied here, namely the analysis of the linearized version of the wave 
equation ([2]) on a fixed Schwarzschild- Tangherlini background. What is more, the 
method of compatible currents being very geometric and robust in nature in fact 
carries over to the non-linear problem suggesting that the decay rates ([6]), (O, dH) 
may be established for the non- linear problem as well. However, in contrast to the 
linear case several non-linear error-terms now enter the various estimates, which 
cannot be controlled a-priori. This requires the introduction of a bootstrap argu- 
ment to be applied in conjunction with the estimates obtained from the method of 
compatible currents. 

It is noteworthy that the paper provides the first application of compatible cur- 
rents techniques in a (non-linear) black hole context. The argument presented here 
is generally more robust than that of [B] but is of course restricted to small data. 

®The quantity l|4]l is related to a boundary term in the vectorfield multiplier estimate associated 
with the vectorfield K. 

^Clearly, this is the associated linear problem to the model of the self-gravitating scalar-field. 
Most notably, it can be treated without any symmetry assumptions on the scalar field, cf. [?]• 
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More precisely, the method presented is expected to be appropriate to eventuaUy 
address non-Unear problems without symmetry, most famously the non-linear sta- 
bility of the Kerr-solution. In particular, since the technique is not bound to any 
dimension one should be able to reprove a version of "Price's law" [Sj for small 
initial data along the lines of the present paper. 

1.3.1 Compatible currents 

The basic idea behind the exploitation of energy currents based on vectorfield mul- 
tipliers is quite simple. We construct a Lagrangian whose field equation generates 
the non-linear wave equation Q satisfied by the squashing field B. The canonical 
energy momentum tensor T^i, can be contracted with a vectorfield to produce 
a one-form P^, = T^,^V^. Finally, Stokes' theorem relates the spacetime (or "bulk") 
integral of the divergence V^Pi, over a certain region to integrals along its boundary. 
This leads to the identity 

/ P^n^ = / V^P^ = / [T^,^^-- + y^V-T^,] . (10) 
Jav Jv JV 

where tt^'' = i [\/^^V'^ + \7'^V^) is the deformation tensor of the vectorfield V . 
One possible application of (jlOp is to estimate a future boundary integral from 
the past boundary and the spacetime-term. On the other hand, for some vector 
fields we will estimate a bulk-term from the boundary terms. The power of the 
method arises from an interplay between the identities associated with different 
vectorfields adapted to the geometry of particular regions. It is crucial that due to 
the Lagrangian structure both the boundary and the bulk term of ((TO)) only depend 
on the 1-jet of B. Suitably applied, the method ultimately produces weighted 
L^-bounds on the fields from which pointwise bounds on the fields follow in the 
standard manner. 

1.3.2 The bootstrap 

Before any bootstrap assumptions can be specified, coordinates have to be defined 
on the Penrose diagram. This turns out to be a rather subtle issue, intimately 
related to the bootstrap argument itself. The crucial observation is that the co- 
ordinates have to be normalized to the future of the bootstrap region, in order to 
capture the decay for late times in the estimates^ This is realized as follows. Con- 
sider the integral curves of the vectorfield Vr, foliating the black hole exterior 
Each of these curves also intersects the curve of fixed area radius r = 2^Mj (with 
M/ being the final Bondi mass the latter is comfortably away from the horizon). 
Hence we can associate a geometric time to any Vr integral curve by using the affine 
parameter along the curve r = 2^JM~f. Now for each such "time" f on the curve, 
we construct a coordinate system Cf (depending on f !) on the black hole exterior 
by the following procedure. We find the point A on the Vr curve associated to f, 
where r^ = 4m. The r = const curve through A will intersect the data at some 
point D. The affine length from ^ to -D along that curve defines the coordinate 
time at A^^ The actual coordinate system (u, v) is finally defined by imposing 
that t = = T holds on the integral curve of the vectorfield Vr starting at the 
point A, at least up to the point B where the integral curve intersects a certain 

^"This is reminiscient of the situation in Christodoulou-Klainerman's proof of the stability of 
Minkowski space |3]. 

^^Notc that for convenience, we have assumed in Theorem 11.11 that the initial data are also 
defined on such a curve, at least up to its intersection with a curve r = rx- 

^^See definition II42I I. We add a factor of \/M in order to avoid dividing by zero when we state 
decay in t. 
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Figure 1: The choice of coordinates. 



constant r^f-curve, fixed once and for all, which is chosen to lie close to the hori- 
zon. Moreover we set r* = '^^^ = at the point A and r„ ~ i (1 — /i) on BE. 
There is some choice to complete the coordinate system by specifying r „ on BC. 
For most practical calculations we will use Eddington-Finkelstein type coordinates, 
setting 1/ = — i (1 — /i) on BC . In any case, the bounds proven will be manifestly 
independent on the choice of u coordinate on BC. 

An important issue immediately arising from the way we define the coordinates 
is that the notion of a constant t slice differs in the different coordinate systems 
depending on the choice of f . (Of course the analogous statement holds for the 
notion of timelike surfaces of constant r*.) Nevertheless, we will show that the 
coordinate systems remain uniformly close to each other in a suitable sense, in 
particular that the coordinate of the initial data slice between tk and the support 
radius is always close to -y/M, however large we choose f. A detailed analysis is 
given in section [8?3l 

Every f defines a T, which in turn defines a region A (T) depicted in Figure El 
It is the region, enclosed by the i = T-curve up to some point B' with coordinates 
(T, rj^.), the null-line v = T ^ r\ linking B' with the horizon, a horizon piece, 
the null-line v ~ + rj^, the t = piece and the u = uq null-line on 

which the field B is identically zero by the assumption of compact support]^. Here 

= supt^2"r* {t,rK). Another curve, r* = r*;, located to the right of r* = 
will also be introduced and fixed. We now choose a small constant c and define 
the bootstrap region to be the region associated to the largest time tb, such that 
for any a/M < t < tb the following "statement "P" holds in the associated region 
A [T (f )) in the coordinate system Cf : 



1. In the subregion {r* > t^} n A{T), the area radius satisfies 



log 



< cVm (11) 



with 



-2V2 - log ■ 



V2 



V2 



and Ma defined to be the Hawking mass at the point (T, r* = 0)0 



(12) 



^^Note that this nuU-hne has a geometric significance by the assumption of compact support. 
The exact value of uq will depend on the coordinate system chosen. 

^^The reader should note that in Schwarzschild with M = the left hand side of l llll l is 
identically zero. The coordinate r* is then the so called Regge- Wheeler tortoise coordinate. 
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2. We haveR 



1 ^ 3 ^ 
-VM < sup t < - VM . 

2 Sn{r*>r*^}n{u>uo} ^ 



(13) 



3. the weighted energy E§ defined in pM)) satisfies (f^ < cM on all arcs 
{t = f <T}r]{r* > r^} C (T). 



4. the energy-flux satisfies m{uhoz,V2) — m {uiioz,vi) < 



for any Ui < U2 



along the part of the horizon located in A (T), where Vi-^- = max (1, Vi). 



m [uri^ ,v) -m {uhoz, v) < 

holds in A{T). Here v+ = max(l,i;). 
6. the integral bound 



(14) 



Fl 



r — — — au < ^ — 



f]2 



for Cl 



sup 



(15) 



holds along lines of constant v in the region {r* < r*i}n{u < T — r* (T, rx)}n 
A (T) , corresponding to a decay of energy as measured by local observers near 
the horizon 1^ 



We define the set 



VT(f) holds in A (T (f )) for all f < f | C 



il/,oo 



(16) 

which will be shown to be open, closed and non-empty. This implies that the 
statement V holds on the entirety of the black hole exterior. The decay rates of 
Theorem 11.11 follow immediately after proving that the coordinate systems used in 
the bootstrap converge to one which is close to the one asserted by Theorem ll.il 

The openness of the set A follows from a straightforward continuity argument. 
The difficult part in closing the bootstrap therefore is to "improve" the statement 
V on the closure of the set A (T). 



1.3.3 Closing the bootstrap 

The third bootstrap assumption is shown to imply -r-i-j- decay of the energy-flux on 

the arcs {t = ti} n {r* > r^-} n {m > jfti}, from which pointwise bounds on the 
field B and its w-derivative are obtained. Additionally, strong decay of B in the 
area radius r can be extracted from the boundedness of . The assumptions also 
provide sufficient control over the coordinate functions at late times. In particular 
one determines the relation between the area radius r (m, v) and the coordinate 
r* = at least in the region where r* > r^. For late times this relation 

converges to the well-known formula expressing the area radius r in terms of the 
tortoise coordinate r* of the five-dimensional Schwarzschild metric as captured by 
bootstrap assumption [T] It follows in particular that the value of r does not change 
much (the corrections are shown to be of order j) along a r* = coTisi-curve in the 

^^This assumption states in particular that the initial data slice is both near and to the past of 
the bootstrap region. It ensures that the bootstrap region does not move away from the data. 

^^That is to say the quantity measures exactly the energy which is not seen by the Hawking 
energy at the horizon. 
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region ?'* > 7'^, allowing us to go back and forth between the two in the course of 
the paper. Moreover, bootstrap assumption [T] is improved. 

Various constant r- and constant r*-curves in the region r > rx will play a cru- 
cial role, since certain integrands arising from the method of compatible currents 
admit good signs in appropriate regions^ The r* = rj;-curve, occurring in the 
bootstrap for instance (along which r w rd by the previous remarks), is determined 
by various requirements defined later but is in any case located to the right of the 
aforementioned r = tk- The latter curve on the other hand, can and will be chosen 
close to the horizon providing a source of smallness in the bootstrap argument. A 
second source of smallness arises from Cauchy stability: After picking some we 
can choose a very late time up to which the fields are still small and after which 
terms like ^lilill^ with C (tk) a constant depending on the choice of tk, are small. 

We now turn to various energy currents arising from vectorfield multipliers and 
describe how the bootstrap is closed. The remarkable properties admitted by the 
Hawking mass for the system under consideration manifest themselves in the iden- 
tity ([To]) for the vectorfield 

The spacetime-term associated to the T-energy identity vanishes and one obtains a 
relation between boundary-terms, which are precisely the associated energy fluxes. 
The monotonicity of the Hawking mass equips all boundary terms with signs when 
applied in the regiorF^ (cf. figure 2]) 



U 



[{h <t<h}n {r* > r* } n{u> uj}) 

{{h +Ki<v<t2+ r*i} n {r* < r*,} n {u < uh}} • (18) 



Such regions arise from a dyadic decomposition of the bootstrap region between to 
and T with <i-|_i = l.lti playing a crucial role later in the argument. 

It can be shown that the boundary-terms associated to the vectorfield 

X = f {r*) [du - d,,) , (19) 

for some carefully chosen bounded function /, are controlled by the energy- flux 
(i.e. the T boundary-terms) and the integral bound ([15]) when applied in the region 
(fT8|) . The function / is in turn chosen such that the spacetime-term of X admits a 
positive sign. In conjunction with the bootstrap assumptions this results in a 

* r* — t ■ 

decay bound for a positive spacetime integral in the dyadic region "'T^^t'^^.'^-^-^, 
which will prove useful in controlling the spacetime integrals of other vectorfields. 

Close to the horizon, in a characteristic rectangle [ui = ti — r*i,U2 = Uhoz] x 
[vi ^ ti + r*i, V2 = t2 + r*i] associated to the dyadic region '^'^"'V'^^^'^J^ , we will ap- 
ply the vectorfield 

Y^^^du+[3{r*)d^ (20) 

for appropriately chosen functions a and (3 (cf. the bold rectangle in Figure[21). The 
strategy is to control the future-null boundary integrals from the past boundary- 



^^By bootstrap assumption [T] constant r and constant r*-curves arc close to one another in that 



region 
18 



From the vectorfield point of view this follows from the fact that T is timelike, that the normal 
to the region is non-spacelikc and the positivity properties of T^^y. Cf. IIIOI I. 
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and the associated spacetime term0 The integrand of the latter contains a part 
admitting a good sign, which can be used in combination with the spacetime term of 
X to control the remaining spacetime term of Y . Moreover, one ingoing boundary- 
term being located completely in the region r* > r*^, is always controlled by the 
energy flux and hence decays like ^. Applying the identity in the characteristic 
rectangle with the bottom being vq = + ^Jz' where an appropriate smallness 
assumption holds by Cauchy stability, and the top being v = v for any vq < v < 
T + r*j^ immediately yields uniform boundedness for both the boundary terms and 
the good spacetime term of Y . The argument can be improved by a pigeonhole 
principle applied in every characteristic rectangle. Namely, one extracts from the 
good spacetime term of F a "good Fy-slice", i.e. a slice on which the local Y- 
energy density decays like times the good spacetime term plus a contribution 
from the ener gy in the region r* > r*^. This is depicted as the dotted line in 
Figure belowr "I Applying the vectorfield identity for Y again in a region with the 
good slice as its past-boundary, one exports the ;jj--decay to all dyadic rectangles. 
Iterating the procedure one obtains decay for all boundary-terms and the good 
spacetime term of Y . The decay of the Y boundary terms leads to the pointwise 
bound < ^ in the region r* < rj;, which can be exported to the region 

^ ^ in* using the energy estimate and the decay in the central region. 
With the pointwise bound on at our disposal, we can finally make use of 

the Morawetz vectorfield 

A = -^du + -^d, (21) 

for a constant As mentioned previously, its application is necessitated by the 
lack of an almost Riemann invariant and it proves crucial in the derivation of decay 
rates away from the horizon. The vectorfield identity for the region =^-'■5^1?!'^'^'' 

associated to any T <T and some large to relates a future boundary term to a past 
boundary term, a horizon-term and their associated spacetime term. 

The boundary terms on the T-arc contain "good" -terms which are precisely the 
strongly weighted energies Eg (t^ of the second bootstrap assumption and error- 
terms. The vectorfield identity is now exploited so as to estimate this "good" term 
on the future arc in terms of all other terms entering the identity. These latter 
quantities are in turn shown to be small or of good sign, which will finally improve 
assumption [3l To derive the smallness for the various terms, it will be necessary 
to subdivide the domain of integration and to apply different estimates in each 
region, carefully taking the geometry of the black hole into account It should 
be emphasized that these estimates belong to the most subtle ones in the paper. 
They make crucial use of the monotonicity manifest in the Raychaudhuri equations 
([23|l and (|24p . and exploit an exponential decay associated with the redshift very 
close to the horizon by introducing an intermediate region between r* = rj^ and 
the horizon. 

For the boundary terms, there are two sources from which the smallness is finally 
obtained: One is the choice of the curve r = r^, which can be chosen very close 
to the horizon. The other stems from the choice of a late time to up to which the 



Physically, the boundary terms of the Y vectorfield correspond to the energy flux as measured 
by a local observer near the horizon. 

^"Alternatively one can extract a "good i^-slice" on which the T-energy flux is improved. This 
will come in handy later. 

■^^This suitably chosen constant defines the origin of the vectorfield. 

■^■^It is here where the pointwise bound on r 2 established earlier enters. 
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initial data has only changed by an amount as small as we may wish by Cauchy 
stability and after which the good decay estimates, i.e. the weight of carries over. 

To establish smallness for the spacetime term appearing in the vector-identity, 
on the other hand, a further argument is needed. This term consists of a "main"- 
term, which is the one that appears in the linear case, and error-terms. The error- 
terms can be dealt with very analogously to the treatment of the error-boundary 
terms. The main term is shown to admit a good sign for r* < and for some 
r* > R* for some R* . The remaining piece in the central region is divided into 
dyadic regions, tj^i = l.ltj. Each i^- integral of such a dyadic region can be 
controlled by tj+i times the spacetime integral of the vectorfield X in that region. 
Since the X-bulk term decays like 



J— — as outlined above, summing up the dyadic 
regions yields smallness for the main fC-spacetime-term (arising from the large time 
to, where we start the dyadic decomposition). This improves bootstrap assumption 

El 

With the third bootstrap assumption being improved on all arcs T < T it follows 
that the decay of the energy has been improved on all arcsl^ As a corollary, the 
same decay is obtained through any achronal hypersurface lying completely in the 
region < r* < -^t. In the final step we find in each dyadic rectangle a "good 




Figure 2: Closing the bootstrap. 

i^T-slice" on which the energy flux is improved to j:;^, very analogous to finding 
a "good i^y-slice" as described above. Combining it with the improved decay on 
the associated arc (cf. the dotted slice in Figure [2|) , the domain of dependence 
property improves the bootstrap assumptions S] and \S\ Additionally, we can finally 
find a good _FV-slice in each characteristic rectangle (improving assumption [5] on 
that slice), which in conjunction with the energy decay now being improved to 
everywhere in r* < j^t, can be exported to all v-slices. Hence assumption [S] is also 
retrieved with a better constant. This completes the proof that the set A is indeed 
closed and the main theorem follows in view of the previous remarks. 

It should be noted that the decay rate that can be extracted in this argument 
is limited by the weights appearing in the K vectorfield, i.e. by the decay in the 
central regionl^ In particular, we cannot derive the stronger decay -^j^ near the 
horizon obtained in [5] for the massless scalar field. It is an interesting question 
whether other methods can improve the decay rates proven in this paper. 

^^This is a consequence of the previously mentioned fact that the expression for Eg contains 
strong weights from which the decay can be extracted. 

Clearly, better decay in the central region could immediately be exported to the horizon by a 
reiteration of the pigeonhole principle in conjunction with the vectorfield Y. 
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1.4 Outline of the paper 



We start by introducing the biaxial Bianchi IX model and some notation (section 
[5]) before defining the aforementioned future-normalized coordinate system Cf in 
section [S] Various a-priori bounds, which can be obtained without invoking the 
main bootstrap argument and turn out to be helpful at many stages of the paper are 
derived in section[3) An important point to keep in mind, however, is that the decay 
of the energy in the area radius cannot be obtained by these methods due to the lack 
of an almost Riemann invariant for the model under consideration. The method 
of compatible currents is explained in more detail in section [5l where moreover 
the relevant identities associated with the regions considered later are derived. In 
particular, the Hawking mass is recovered as a potential of a certain vectorficld- 
currcnt (section [6|) . After defining the bootstrap assumptions (section [T]) , various 
bounds for the fields are derived from them and the stability of the coordinate 
systems Cf defined in section [3] is established (section [8]). The identities associated 
to the vectorfields Y and X are analyzed in sections |9] and [TOl Here a somewhat 
lengthy argument is pursued to construct the function / implicit in the vectorfield 
X, which finally ensures that its spacetime term admits a positive sign. Section 
111! reveals how to control the weighted energies produced by Y near the horizon 
with the help of the vectorfield X. The relevant version of the pigeonhole principle 
is also explained at this stage. Finally, in section [12] the Morawetz vectorfield K 
is introduced and the necessary estimates to control the various error-integrals, as 
outlined in the introduction, arc performed. Everything is put together in section 
1131 where the bootstrap is closed. The paper finishes with some final remarks and 
open questions. 



2 Biaxial Bianchi IX 

The class of biaxial Bianchi IX metrics was introduced in [1] . We recall that these 
spacetimes are topologically M. = Qx SU{2), where Q is a two-dimensional mani- 
fold and that global coordinates (it, v) can be found on Q expressing the metric of 
M in the form 

g = -n^ (u, v) dudv + {u, v) (e^Biu^^) (a? + a^) + e-4B(«..)^2^ (22) 

where B and r are functions Q ^ M and the cr; form a basis of left invari- 
ant one-forms on SU{2). Note that if i? = 0, the symmetry is enhanced to 
(SU (2)^ X SU {2)j^) /I? = SO (4) and the metric reduces to the five dimensional 
Schwarzschild-Tangherlini metric in view of the higher dimensional version of a 
well-known theorem due to BirkhoffEl In this sense, B is the dynamical degree of 
freedom ruling the model. See [5] for a more detailed discussion. 

The vacuum Einstein equations for the above model reduce to a system of 1 -|- 1 
dimensional PDEs on the quotient manifold Q: 

a„(l]-2a„r) =-^((5,0'), (23) 
9„ {^-^d,r) = - J {{B.^Y) , (24) 

6 r r r-^ S r \ 2 J 



^^Note also the relation between the familiar round metric (dui^^) and the bi-invariant metric 



on duj^^ = i (aj + a^+aj). 
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dud. logo = g + Ir,.. . - 3 (B,.) (B,^) = + |i (p - - 3^ , (26) 

B.u. = -I^S,. - + S (e-«^ - e-^^) . (27) 



Here we have defined the quant it-J^ 

P = 2e-^« - ie-«« < ^ , (28) 

with the inequality following from elementary calculus. Equality holds if and only 
if i? = 0. Note that the non-linear wave equation (j27p can be written as ^ with 
□ being the d'Alcmbcrtian of the metric ([22|l . 

A remarkable feature of the above system is the existence of a function m (w, v) 
called the Hawking mass and defined by 

Since the inequalities < and > were shown [5] to hold everywhere on the 
black hole exterioj^, the Hawking mass has the following monotonicity properties 
there: 

1--, 



dum = -4r3— (B „)' +riy[l--p] <0, (30) 



d,m = -4r3^ {B,„f + rX (^1 - ^p^ > . (31) 

This allows the derivation of energy estimates for the field B, which plays an im- 
portant role at all stages of the present paper. The existence of these estimates was 
already an essential ingredient in the proof of the orbital stability [S] . 

We conclude this section recalling some notation introduced in . We set 

A = r,„ j/==r„ C-^-^S.^ d^riB^^ (32) 
and introduce the quantities 

and 7 = = — (33) 



I - p -4i/ 1 - M 4A 

satisfying 

2C' 



(34) 



2 6^ ~ 

^,v=l[ 3--- ) , (35) 



r2 A 



as well as the auxiliary quantities 



^i{B)=[l- 8B' and ^,{B) = ^ - e''^) + 8B^ , (36) 

both of order B^. The volume element associated to ([^ is 

dV^oZ = y/gdudvdw = —dudydA^s = r^Vl^dtdr*dAs3 (37) 



The quantity p is related to the scalar curvature of the group orbit by R= ^ p. 
^''They hold on the initial data for small perturbations of Schwarzschild-Tangherlini and are 



seen to be preserved by equations iTAl and I l24t . 
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where in the standard Eulcr-coordinates on SU{2) (cf. [5]) 

11 f 

dAgs = —sindduj ~ — sm 9 dOdcftdip and hence / dAgs = 2t:^ . (38) 
8 8 J 

The monotonicity of the Hawkmg mass justifies the definitions 

fTT-mim "fTT-max for the minimal and maximal Hawking mass. (39) 

Furthermore, the quantity Mf will denote the final Bondi mass and M the mass of 
the perturbed Schwarzschild solution. The mass M determines the scaling of the 
problem and I have normalized all quantities appropriately using factors of M. In 
particular, "smallness" always refers to dimensionless quantities. 

We write C (e) for a constant satisfying limg^o C (e) = 0. The notation a ~ & is 
used if there exist uniform constants ci, C2 with ci < ^ < C2. Finally, we define 

u+ =max(l,w) and Vi+ =max(l,?Ji) . (40) 



3 Choice of coordinates 

As mentioned in the introduction, the choice of coordinates is already a rather 
delicate issue for the problem under consideration. Although the final result does 
not depend on the choice of coordinates, the bootstrap-techniques applied in the 
proof require the coordinates to be normalized to the future of the bootstrap region 
introduced in section [71 If on the contrary one normalized the coordinates on the 
initial data, one would not be able to obtain the improved decay of the fields at late 
times from the estimates, roughly speaking because contributions from the initial 
data, which have not yet decayed, enter the estimates. This necessitates, after 
a purely geometric definition of "time" for V?' integral curves on the black hole 
exterior, the introduction of a different coordinate system Cf ~ {uf,Vf) defined 
with respect to every such "time" f . All such coordinate systems Cf are defined on 
the set 

v = j- (1+) n J+ (^Sr^) (41) 

of the black hole exterior. In section 9 we shall exploit the bootstrap assumptions 
to establish that - in a certain region - these coordinate systems arc uniformly close 
to each other in a suitable sense. It should be observed that the coordinate systems 
Cf are different from the coordinate system asserted in Theorem 11.11 In the last 
section of the paper we will show that the coordinate system Cf for f — > oo is close 
to the one asserted by Theorem ll.il 

We begin by considering the family of Vr integral curves starting out from some 
r = rK-cmve which is chosen close to the horizoro such that still 1— /i > c > holds 
for a small c, and ending at spacelike infinity i^. These curves foliate I?n {r > r^}. 
Moreover, every curve admits a unique point where r = 2-\/m. We pick any such 
curve and label the corresponding point by A. Denote the mass at A by tua and 
consider the curve = ArriA going through A and intersecting the initial data at 
some point D. Let tad be the affine length of the constant r curve (with tangent 
vector normalized to one) connecting A and D. Finally, define 

T^Vm+ '^^^ =Vm + V2tad (42) 



to be the time associated to the Vr curve under consideration. In this way we 
can assign a notion of time to any Vr integral curve. Considering next the curve 



*The choice will provide a source of smallness later in the bootstrap argument. 
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4Af f with affinc parameter f starting from the initial data, we obtain a map 
'd:[0,oo)3T^Te\VM,oo) (43) 



which is defined by taking f to the time associated to the Vr integral curve which 
intersects the curve = AMf at f. The map d is easily seen to be continuous and 
surj active. 

For every f a coordinate system (m, v) is defined as follows. Let A have coordi- 
nates (m, v) — [T = § [f) ,T = § (f)). Set K = 7 = ^ along the Vr integral curve 
up to r = rx- Since 

Vr {u + v)^ {Vrr + {Vrf = A „ A) = ^ii^ (7 - , (44) 

we have that t = (thus defining t) is indeed equal to the constant T on the 
Vr integral curve through A. Moreover r* = ^^^-^ is equal to at A. Let the Vr 






Schwarzschild 



Figure 3: The choice of coordinates. 

curve defining the coordinate system intersect r = r/^ at B. We erect the constant 
u = tiB-ray to the past of B and set k = ^ there. The coordinate system is 
completed by specifying the u-coordinatc on BC. We set = — i (1 — /i) on BC. 
This might send the horizon to u = 00, namely if 1 — /j, = at Ccj We will see that 
in these coordinates w ^ 00 at The coordinates thus defined will be refered 
to as Eddington Finkelstein coordinates. We also allow ourselves to move freely 
between the coordinates {u,v) and [t = ^^^,r* = ^^^)- 

Clearly if f = 0, then the associated integral curve coincides for r > rx with 
the curve on which the initial data is defined, and t ~ \[M defines the initial-data 
slice in r > tk- Note that in any coordinate system associated to some f > 0, a 
slice on which t = constant does in general not agree with a Vr-slice. However, 
once we have introduced the bootstrap assumptions, we will be able to show that 
the two slices mentioned remain uniformly close to each other in r* > rj^ for any 
f > 0. This argument is postponed to section [8.3.21 Here we only introduce 

Notation 3.1. Let t^^ denote the t- coordinate, measured in the coordinate system 
defined by ta, of the point defined by the intersection of the Vr integral curve 
determined by tb and the curve r^ ~ Am (§ {ta) , r* =0). 

We conclude with a remark on the differentiability of the coordinate systems. 
Due to the "cusp" at the point B the coordinate system is only : The quantities 
K and 7 (and by definition ([33]) the first derivatives of the area radius function 



^^Of course, one does not expect this to be the case gcncrically. 
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r (u, v)) are clearly continuous. The second derivatives r ,j„ and r^uu however, are 
discontinuous at the point B. This could be avoided by applying an appropriate 
smooth interpolating function in a small neighborhood around the point B. How- 
ever, we will see later that the bootstrap involves only first derivatives (and hence 
continuous quantities) and that the regularity suffices to close the bootstrap. 



4 Basic estimates 

In this section we are going to show that given an appropriate smallness assumption 
on the field B, namely ^ and the field and its derivatives remain small on the 
entire P. Since this "first round" is independent of the main bootstrap argument, 
it provides a good way-in to familiarize oneself with the basic estimates applied 
in different regions of the black hole exterior. The bounds in this section will be 
proven in the coordinate system Cf associated to any f > 00 In this context it 
is crucial that the smallness assumptions ^ and ([5]) are manifestly independent of 
the coordinate choice. From [5] we recall that 

777 

1 - — ^ < e {5) with lim e ((5) = (45) 

nrimax -5^0 

can be chosen arbitrarily small by an appropriate assumption on the initial data. 
We will abbreviate e [5) by e in the following. In view of the monotonicity-properties 
of the Hawking mass ( ([5D|) and pTjl ) the mass difference between any two points 
cannot exceed m„iax • e (^)- We note 

Lemma 4.1. // ^5\ l holds, then on the horizon we have 

< 1 - M < (46) 

Proof. From [5] we have both 1 — > on the black hole exterior, as well as the 

Penrose inequality 1 J- < holding on the horizon with Mf the final Bondi 

mass. Combining this with ((45|) immediately yields ([46|) . □ 

Corollary 4.1. The area radius r satisfies 

|^^_^_|<ip^e onH+ (47) 

with r± being the maximal (minimal) value of r on the horizon. 

Corollary 4.2. For any given rj > we can choose the S of the initial data so small 
that for some r ~ rx curve located completely in V the estimate 

rx - r < rj (48) 

holds in r < rx . 

For the estimates in this section only we will explicitly couple the location of 
the r = rx curve to the smallness of the initial data. In particular we define the 
curve rx by 

1 - = e3 (49) 

It follows easily that the maximum r difference in the region r < rx satisfies 

Ar<rx-r.< S-^^ei . (50) 



^^Note that if r = the coordinates are normahzcd on initial data. 
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Proposition 4.1. In any coordinate system Cf and with the assumptions of Theo- 
rem \1.1\ on the initial data we have 

r\B\ + ^\e\ + M-^ ^ < \/m ■ C {5) (51) 

V 

everywhere in T) . Moreover the coordinate junction n satisfies 

1 



2 



<C{5) (52) 



everywhere in T) . Here the constant C {5) can he made arbitrarily small by choosing 
the 5 of the initial data sufficiently small. 

Proposition 14.11 will immediately follow from Propositions 14.2114.51 (plus their 
associated Corollaries) proven in the remainder of the section, each of them estab- 
lishing the bounds in different regions of the black hole exterior. Note that the 
radial decay of B promised by Proposition 14. II is weaker than that of Theorem ll.il 

Proposition 4.2. In the region T) we have 

\B{u,v)\<Ci{e{6),5) . (53) 

where Ci (e (S) , S) can be made arbitrarily small by choosing the 6 in (0) small 
enough. 

Proof. We integrate from the initial data to any point in the region r > r^ and 
estimate as follows 

, 3 



\B{u,v)\ < si iudata,v)+ / —du 



Udata 



'\/m\ ' 1 



< S + -^\Jmmax - VTlmin SUp . , 

\ r V2 r>rK Vvl - A*/ r 



which proves (|53p in that region. 

3 

Next we turn to the region T> r\ {r < rfc}. We choose a constant C > 2M^r_^S 
such that \B\ < C still implies that 

^-p=^-(2e-^-^-^-e~^A<^^ (55) 
2^2^ 2 y-2r2 ^ ' 

in 2?n{r < rx}- For tk sufficiently close to the horizon it is easily seen that C ~ 
is good enough. Define the region 

i^{u,v) eVn{r <rK} : |-B (m, v) | < C for aU (m, w) e J" (u, w) | (56) 

which is clearly open and non-empty. We are going to apply a bootstrap argument: 
Pick a point in the closure of TZ, where B < C hy continuity. We are going to 
improve this bound by showing that in the causal past of that point B is in fact 
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smaller than ^. By continuity it follows that the set TZ is also closed, hence must 
constitute the entirety of I?n {r < rx}- The argument proceeds in two steps. First 
we make use of the redshift estimate integrating the equation 



4 K 



2r 3 

from the initial data yielding 



C /4k 



4k 
37 



and hence 

c 



36* 4 k 



-8S „-2B\ 



2r 3 



(57) 

(58) 
{u, v) dv 



(w, v) 



< 



g -ffii'^™]^"''"'''^'" K {u,v) dv\l I — ^{u,v)dv 



4 

o sup 



. -8i3_g-2S| 



1 



g is [7S'™]'^"'^)''^K (u, w) di; 



< 



3 ./rZuir, 



(59) 



It follows that 1^1 is bounded (but note that the last term might not be small) in 
that region. In the second step we integrate from the r = tk curve, on which B is 
small by (|54|) . or the initial data to obtain 



(60) 



r ( - 

B {u,v) ~ B {ujf ,v) = / —{u,v)du 
and use the previous bound ([54]) 



\B{u,v)\ < S 3— + £3, 



< 6- 



C Mi 

2"(Zy 



3" 

7-2 



y((rK)^-(r_)^) 



Now because the r difference is given by (|5G|) in the region under consideration, we 

c 

2 



have indeed shown that B is smaller than ^ in 7?. for an appropriate choice of e. 



By continuity the set TZ is also closed. Hence TZ = V H {r < tk}- 
Corollary 4.3. In r > rK we have that 



\B{u,v) I < VM- 



Proof. This is the statement of ([5 



□ 



(61) 



□ 



It is instructive to compare the ^-decay of Corollary 14.31 with the analogous 
estimate derived for the massless scalar field in four dimensions [6]. In the latter 
case, one obtained by the above method -^-decay. There existed an almost Rie- 
mann invariant, i.e. a certain combination of the field and its derivatives, however. 
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admitting better decay properties than the field or its derivatives alone. Via this 
quantity, it was possible to improve the decay in r of the field itself to ^ , which was 
it turn sufficient to extract energy decay in r. In five dimensions there is no almost 
Riemann invariant and energy decay in r will only be obtained from the application 
of the Morawetz vectorficld K in the context of the bootstrap argument pursued 
later. 

Corollary 4.4. In the region TZ — T) f] {r < tk} we have 

^ <M-^C3ie,S) . (62) 



Proof. This follows from revisiting the red-shift estimate (|59[) above, this time im- 
proving the estimate for the (e^^^ — e^^^)- term by Proposition |421 which implies 
that \e~^^ — e~^^| is e-small. In this way we obtain a smallness factor for all the 
terms involved in (15911. □ 



Proposition 4.3. In T) we have 

^ ' < C4 (e, 5) . (63) 



Proof. Integrating p4p from the t = Tr\{r > r^} surface to any point in the region 
{r > rji} yields 

/ 2 (1 - /i) \ _ 
k{u,v) = K{uT,v)ey.Y>\ I ^-^-q -]{u,v)du. (64) 

If the point under consideration lies to the future of the t = T hypersurface {u > ut), 
the upper bound n < \ follows from monotonicity, whereas the lower bound is 
obtained via 

, / M 1 ( ( 2 \ r c' (1 - /^) . ^ , 

\k.{u,v)\ > -exp sup ^— / {u,v)du 



1 / — 4mn 



> ^cxp -=-) . (65) 



On the other hand, integrating (|34p from the null line u = T — r* (T, tk) downwards 
the lower bound follows from monotonicity and the upper one by using 



|«;(u,«)|<iexp((C3(e,5))')cxp(-^^ + ^) <\ + ~c^{e,5) . (66) 
Z \ / \ r yu, V) rK I 2 

Since the r-difference in the region r < rx is e "J -small by (|50p . we obtain the desired 
upper bound for k in particular on all of r = rx- 

Now any point located in the past of the t = T hypersurface and satisfying 
r > rK can be reached by integrating ([34]) from either t ^ T or from r = rK where 
the upper bound ([55]) has already been established. The lower bound for k at such 
a point follows from monotonicity, the upper one from 

|. (., .) I < Q + C4 (e, S)^ exp (^sup £ v) du 

< (i + 5. (.,.)) exp [il^^). (67) 
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To extend the estimates to the entire region r < tk '^c integrate (j34p from the 
r = tk curve (on whieh the lower bound (|66p and the upper bound (|65p has been 
estabUshed) to the horizon. Again the upper bound follows from monotonicity and 
for the lower one we write 

K {u,v) — K {uq,v) cxp I / — :^\{u,v)du (68) 

\Juo (l^) J 

and estimate, using ([5^ 

\k{u,v)\>^ exp I I exp ({C3 (e, d)f^ 



exp 



(69) 

Taking again ()50|) into account, we obtain the lower bound for k also in that region. 

□ 

With the bound on k established we also have good control over the quantity 
A = K (1 — /i). In particular A < 1 everywhere and A becomes very small (perhaps 
zero) at the horizon. In particular, it follows that 

holds everywhere on >§n {r > rx} and hence the -l-part of the smallness condition 
([3|) implies smallness for as well. With this in mind we can prove 



Proposition 4.4. In T) we have 



\0\<C,{e,6)J- (71) 
V r 



Proof. We rewrite equation (p7|) as 

and integrate it from 5 n {r > rx} to any point in V. We note that for \B\ small 
we can find a constant K such that 

{e-'^^e-'^)'<K(l-^-p^ (73) 
holds. This constant approaches | as \B\ goes to zero. We then estimate 



\0[u,v)\ < Mi--(Udata,v) + — du. — — (u, V) du 



+ SUP ^-.j ^K^j^^ji-.) [1 - -pj - in,v)dn 

<Mi^ + ^^^+8Vi^.(74) 

□ 

Finally, we extend the bound on ^ to the region r > tk- 
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Proposition 4.5. We have 



in all ofD. 



<Ceie,S) 



(75) 



Proof. Integrate equation (j57p from the r = r/<--curve, where |-| < C3(e,(5) by 
CoroUary 14.41 out to infinity. Note that due to tlie estimate proven for the field B 
in Corollary 14.31 we may achieve (choosing S small enough) that 



3 3m 



(76) 



2 

holds in the region r > rx- Using again (j73p wc can follow the string of estimates 



- [u,v) 


< 


- (M,WrA-) 




'L 


V 




V 





< 



+ 



< 




to conclude the result 



So far we have shown that rB, ^, and y/rO are small and that k is everywhere 
close to ^ for the perturbed spacetimc. Estimates for some higher derivative quanti- 
ties will be required later. However, since all bounds can be considerably improved 
once the bootstrap assumptions have been introduced, wc postpone the derivation 
of further pointwise estimates to section [8.41 Here we only note 



Proposition 4.6. On V we have, independent of the coordinate system Cf, the 
bound 



m 



< 



1 



Crie) 



(78) 



Proof From the fact that k = ^ on {t = T} n {r* > r* {T^rx)} (hence n.^ 

there) and on {u = T — r* (T, r^)} D {t < T} (hence k „ = on this null-line) the 

n , 







bound ([78]) follows on these sets. Wc can obtain the quantity at any point on 
T) by integrating equation (|26p from the aforementioned set to the desired point. 
Inserting the estimates of Proposition 14. II gives ((78)) everywhere. □ 



Remark: The quantity is discontinuous at the point B in the coordinate 
system Cf. This discontinuity is propagated along the null- line v ^ v (B) when inte- 
grating the quantity du-^ (which is continuous! (cf. I26[) ) in u (cf. also Appendix lX)) . 

We conclude the section with a useful bound for the quantity 7 in the region 
Vr\{t<T}C^{r> rx}- 
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Proposition 4.7. InV n {t < T} H {r > rx} we have in the coordinate system Cf 



C8(e,<5) <7-^ <0 



(79) 



Proof. Integrate ([55)) from the t = T-slice in the past-direction. By monotonicity 
7 < ^ is obvious. The other direction is derived from 



7 (u, v) — J {u, vt) exp 



(80) 



and the estimate 



7 (w, v) > - exp 



sup 



{r>rK}n{t<T} (1 - M) 



Vt 02 

— (u, v) dv 



> -exp[C8(6)], 
which follows by choosing the mass fluctuation small enough. 



(81) 
□ 



We close the section by emphasizing once more that the bounds proven in this 
section are independent of the particular coordinate system used, i.e. of how large 
we choose r (and hence T). In this context it is important that the smallness 
assumptions ([3]) and (O are invariant under a change of coordinates. 



5 Compatible currents 

5.1 The basic identity 

Varying the Lagrangian 

>C = Ig'-'^d.Bd^B + ^ (l - ^p) (82) 

with respect to B leads to the non- linear wave equation ^ satisfied by the field B. 
We associate to the energy momentum tensor 

T^, = d,,Bd,B - ig^., {dBf - ^g^, (^1 - (83) 

satisfying the equation 

V^T^, = l(^l-|p^V,r. (84) 
Given any vectorfield V we can define its deformation tensor 

= I (^"V" + VV^") (85) 

and the vector 

pa ^ g°'f^TfisV^ . (86) 

The method of compatible currents is based on the following basic identity for an 
arbitrary vector field V: 

- V„P" = - (Tapnf + (V^T,^) l^") . (87) 
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5.2 Useful formulae 

In (u, u)-coordinates the components of the energy momentum tensor (|83p read 



T 


= {duBf , 


T 






1 


T 






1 



2 \ 1 / 2 



The vectorfields used in this paper have u and v components only and will 
furthermore depend only on these two variables. For such vectorfields we compute 
the components of their deformation tensor: 



= -'feK + AK). (89) 



Finally, the following explicit formulae for the contraction 

T^.TT^'' = T„„^"" + T™^''^' + 2T„„7r"" + TyTr''^' (90) 

will be useful: 

'^^K + Ak) + (91) 



r 



and 



5.3 Basic regions 

In the course of the paper we shall apply the basic vectorfield identity ([57|) for 
different vector fields in adapted regions of the black hole exterior. Here the relevant 
formulae arising from ([57)) for these regions are derived!^ 



^^Since the coordinate system is only picccwisc , the justification of those formulae, which 
are easily derived formally, requires some care. A detailed discussion can be found in Appendix 

El 
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5.3.1 Characteristic Rectangles 

Writing out the identity ([57)1 for a null- rectangle TZ ~ [ui,U2] x [wi,W2] yields 



vol 



Defining the bnlk term 

ivoi{n) 

and the boundary terms 

?v 



/ / / [du{^P'') + d.,{^P'^')]dudvdoj. (93) 



[Tc^pTrf + (yf'T^p) 1/" j —r'dudvdAs^ 



F^i[u,,U2]x{v}) = - 



y/gP^ (u, v) duduj 



r-" {duBY V + — [1- -p]^ 



(94) 



dM,(95) 



Fl{{u}x[v,M) 



we find the identity 



y/gP^ {u, v) dvduj 



dv , (96) 



F^{{U2} X [vi,V2]) + F^ {[UI,U2] X {V2}) 

II (TZ) + F^ {{ui} X [vuv2]) + Fl {[ui,U2] x {^;i}) 



(97) 



5.3.2 The region """V^^^/^-^ 
Another important region is 

""^^'".1 < ^ < ^2} n {r* > r*} n{uj<u< uh} 

U ({(u, v) e [ti - r*, uh] X [ti + r;, t2 + r*] }) 
for which one finds the basic identity 



fv = it2) - F^ (t,) + , 



with the bulk term 



and the boundary terms 
1 



[tl,t2l 



(-r^,^(^'^-(v^r^,)y")dyoZ 



27r2 



i^^ (0 = 



-P* {t,r*)n'^r^dr* 



Uh 
t-r* 



r- {d^BrV- + ^-^(l-lp]V^ 



(98) 



(99) 



(100) 



du, (101) 
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where 



yu 



2 (duB) 



2r2 



yv 



2 (d^B) 



2r2 



and 



and 



1 



27r2 



1 



J'' 



27r2 " 



2t2—Ui 



2ti-uj 



(102) 



(103) 



(uj,^)^^. (104) 



For the region under consideration we will also need to apply Green's identity to a 
term of the form D ■ □ (-8^) for some function 



+ j [{nB^)D]dVol^ ... + j [B^{nD)\dVol 

+ G{t2)-G{h) + N{t2)~N{h) + H^^-J^,{im) 



where 



1 

2^ 



[B'^dtD - DdtB'^] {t, r*) dr* , 
^A^(<)= J [B^duD- DduB^]r^u,t + r*)du, 
I [B'^dyD - DdyB^] {uh, v) dv , 

Jti+r* 



27r2 
1 



27r2 
27r2 



[B'^d^D - DdyB'^] {u,j, v) dv . 



(106) 
(107) 
(108) 
(109) 



^■^The formula derived here is a- priori valid only for D £ . However it also holds for a D 
admitting less regularity, as is shown explicitly in Appendix [A] whore we demonstrate that for 
the cases where HIOSI I is applied in the paper (equations (137311 and II260II ). D indeed satisfies these 
requirements. 
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Wc then define the rcnormalizcd bulk term 



■= ■ + i .J.. (D^)] dVol , (110) 



[tl,t2l 



for which the identity 



with 



= ih) - (to + H.,, J„, (111) 

F^it)^FUt)^Git)^N{t) , (112) 

Huh = ^^ff - Hufj , (113) 

= Jl - Juj (114) 

holds. Note that for uj = uq, the boundary terms Juj all vanish, because i? does 
not have any support on u = ug by the domain of dependence property. 
Finally, for future reference we also define the subregion 

and the slice 

= f{< = t}n {r* > r*i}) U ({v = f + n {r* < r*,}) . (116) 



6 The vectorfield T and the Hawking mass 



Recall that the Hawking mass m defined in satisfies ([50]) and (PT|) . The one- 
form dm is closed and by simple connectedness of the Penrose diagram, exact. It 
follows that energy is conserved. This fact can also be seen from the integral identity 
([87)1 applied to the the vector-field 



_ 4A 4i/ 



(117) 



If we apply the identity ([871) in the region ""^^jj'^'j^j, energy conservation translates 
into the following relation between the boundary terms: 



Fb (t2) = FE (ti) - H^^ + J^iuj) , 



(118) 



where 



Flit) 



t-r 



cl 



(it, t + r*i) du 



4^ V (^-"^ +r{\~v)\\--p] I (t, r'^) dr* , (119) 



{uH,v)dv, 



(120) 



2t2 — 'tij 



(uj,z;) . 



(121) 



We will sometimes use the notation E (E), for the energy flux through an achronal 
slice E. 
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7 The bootstrap 



The bootstrap is intimately related to the choice of coordinate systems defined in 
section [3] We will use the notation introduced in that section. 

7.1 The bootstrap region and the statement V 



Let 



-3\/2 - log 



'2 - 



and c be some small constant. Define 

S^tdt + (r* - a) dr" 

and the quantity 

27r2 



2 + ^/2^ 
S = tdr* + (r* - a) dt 



M 



with 



+ (-2.)( [SB + l'—^B 



t<T 



) {(SBf + (SB)'- 



:B' 



(122) 



(123) 



(124) 



To each r we associate the region A (T (f)) = 

V V // [2VM,T) 



(125) 

(hence defining the T 




Figure 5: The bootstrap region. 

in ^25]) ). 
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Wc define the statement 7^t(t) associated to a region A {T (f)) to btj 
1. In the subregion {r* > r^} n A{T), the area radius satisfies 



with 



r(i,r*) + \/^(log 



P 



-2V2 - log 



2 - v/2 



< cVM (126) 



(127) 



2 + 72 

and Mji defined to be the Hawking mass at the point (T, r* = 0). 

2. We have 

i\/M< sup t<|VM. (128) 

3. The weighted energy-density (|124p satisfies 

1 



M 



4. The energy-flux satisfies 



Ej^ [T] <c on aU arcs {2VM < t = T < T} n {r* > r^^} n A (T) 



M 



m{uhoz,V2) - m{uhoz,Vi) < c M- 

(wi+) 

for any vi < vi along the part of the horizon located in A (T) and 



(129) 



(130) 



m (ur*^ , d) - m (w/ioz, v) <c 

holds in A (T) for an r*j defined in the subsection below. 
6. The integral bound 



r^^^^^du <ClM^ for Cl = sup 

-1/ r*>r*^ 1 - M 



(131) 



(132) 



holds along lines of constant v in the region {r* < r*;}n{it < T—r* (T, r^^)} C 
^ (T), corresponding to a decay of energy for local observers near the horizon. 



Finally, we define the set 



A 



|f e \/M,oo^ VT{f) holds in ^(T(f)) for aU f < f I C ^/M.oo 



(133) 

Note that the lower bound on f ensures that T > 2\[M (cf. (|42p ). The following 
key- Theorem will close the bootstrap and is easily seen to imply the decay rates of 
Theorem ll.il It will only be proven at the end of the paper. 



Theorem 7.1. The set A is non-empty, open and closed. 



■^^We will sometimes abbreviate T (f ) by T, reminding the reader that any T arises from f as 
described in section |3] 
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A few remarks are in order. The first two bootstrap assumptions ensure that the 
diff'erent coordinate systems Cf do not move too far away from one another, at least 
in the region r* > r^. The first controls the deviation of the relation between the 
coordinate r* and the area radius r from the familiar relation between the Regge- 
Wheeler coordinate and the area radius in the Schwarzschild metric. In particular, 
for Schwarzschild the left hand side of (|126p is zero. The second assumption ensures 
that the bottom of the bootstrap region (the t ~ slice) does not move away too 

much from the geometrically defined initial data (and is moreover always located to 
the future of the data). In other words, the coordinates of S'n{r* > fjf }n{u > uq} 
are similar in all coordinate systems Cf. 

The open-part of Theorem 1 7 . 1 1 follows from a simple continuity argument: 

Proposition 7.1. The set A defined in US 3]) is open. 

Proof. We observe that the integral Eg (t) and in fact all the quantities appearing 
in statement V of the bootstrap assumptions depend continuously on the choice of 
f. ' □ 

One should note in this context that all bootstrap assumptions involve only first 
derivatives of the fields and the area radius, and hence only continuous quantities 
(cf. tlie remarks on the differentiability of the coordinate systems at the end of 
section [3|) . 

The hard part of Theorem 17.11 consists in showing that A is closed. This will 
be accomplished by improving the constants appearing in the inequalities of the 
bootstrap assumptions. 



7.2 The choice of r*i 

In this subsection we define the quantity r*i with respect to the coordinate system 
associated to the bootstrap region. Clearly, the location of r*^ will change between 
different coordinate systems when the bootstrap region is altered. However, by 
bootstrap assumption [1] it will always stay close to a geometrically defined curve 
of constant r, which is determined below. 

By Propositions 14.61 and 14.11 we know that on V the bound 



^ .1' 



m 



1 



<Cie) 



(134) 



holds in any coordinate system Cf. For any small number rp > we can hence 



choose the initial data small enough such that there exists an ry < M satisfying 



max 

r<rY 



ry 0„ 1 



(135) 



Here CoroUarv 14.21 has been used for the bound on the first factor. By bootstrap 
assumption [1] the curve ry ■= infto-r* (t^ry) is always close to the geometrically 
defined curve ry. Hence we can additionally impose that 



M 



(136) 



holds. Next we are going to determine how small i]j has to be. We define two 
hmctions a (r*) and (3 (r*) in the coordinate system associated with the bootstrap 
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region as follows. 





The function a which is supported only for r* < —^^/M is everywhere non- negative 
and defined by setting a ^^r^ = t+i j^--u,,a^ ^ _ 



r 



a (r*) 



/M 



for r* < 
in [r^,r*j^] 



M3 



(137) 



^ (\/M+|r*|)2 



T in[r^,ryj 



with M = m (T, r* = 0) and x a smooth positive interpolating function. In partic- 
ular a = 1 on DC. 

The non-negative function /3, again with support only for r* < — i V^M, is defined 



by setting (3 [r}j = 



24 



= and imposing that 
18 



-n^it,rn>P'>——n^it,r*) 



r{t,r*) r{t,r''^ 
in all of ?'* < ry ■ We can estimate the value of (3 on by 



(138) 



P^r*dr* < I 24 — dr* < 



2A^^dr* logrdr* < 12 log — 

7 -f K 



Hence /3 remains controlled by the r-fluctuation in r* < ry and hence small by 
choosing tjj above suitably small. Note that a and /3 are in particular supported 
away from the curve r* = 0. 

We finally choose the ip of (|135p . (|136p so small that the inequalities 



4a— — a'r ] > max 



2 I —a-(3\ 
4k 



n. 



r K (1 



a > K I 4/3A + 2r/3' -t- 8r/3-^ ) -I- max 



4VM 



/A/ 



2VM 2VM 

> 45 



24^r-^ -f (1 - /i) (-70k - 36k/x) 



hold in the region r* < ry and set r*j = — 2\fM . 



(139) 

(140) 
(141) 



Remark: The constant ?/; and the corresponding ry (and the upper bound on 
initial data) can easily be computed explicitly and is fixed once and for all. In 
particular it does not depend on the size of the bootstrap region and the coordinate 
system that comes along with it. The curve r* = ry and hence r* — r*^ is then 
also fixed and always close to ry by bootstrap assumption [1] and the fact that tk 
is chosen much closer to the horizon than ry. Smallness for the bootstrap on the 
other hand, will be exploited via the r^-curve and by choosing the initial data even 
smaller to "beat the constants" which are introduced by the choice of r*j . 
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7.3 Cauchy stability 

For the closed-part we will have to improve the constant c in the statement V 
(i.e. the bounds (|129m32|) ) in the region A (T). The argument constitutes the body 
of the paper. In this context we note that within the process of improving the 
bootstrap assumptions there will be two sources of smallness. The first arises from 
the fact that r = rj^ can be chosen very close to the horizon. The second is obtained 
by selecting a Vr-slice belonging to some large tq (and hence large associated time 
to) up to which Cauchy stability holds by a suitable smallness assumption on the 
data. This is expressed precisely by the following 

Proposition 7.2. For any small rj > 0, 5 > Q, and any large tq (hence large 
associated time Tq = (tq), with 1} defined in we can find an rn and a d > 

such that the following statement is true: If the smallness assumptions and (0) 
of Theorem ] 1.1\ hold for S, then 

1. The curve r ~ r^ away from the horizon satisfies rj^ — < 77 

2. in the coordinate system defined by f, G [0,fo] the t-coordinate of the subset 
S n {r > rji} n {li > uq} of the initial data satisfies 

\t - VM| < <5\/M. (142) 

3. In the coordinate system defined by Cfg , the statement V holds with constant 5 
(instead of c) in the region ^'""''^^^2^1^ t ] ^"''^ moreover, the pointwise bound 



\B\ + 



holds on any slice St (cf. U16\) ) for 2yAI < < < To. 



M~i\e\<y^ — (143) 



Proof. The first assertion is the statement of CoroUarv 14.21 For the second state- 
ment consider the coordinate system Cr^fl(f^) for a given f, € [0,fo]. The vectorfield 
Vr introduced in section [3] can be expressed in the associated {t,r*) coordinates 

'^r=-^[{-f~K)dt + {-/ + K)dr*] (144) 

4^7 

as can the vectorfield \/±r which is defined to be orthogonal to Vr and whose 
integral curves are the curves of constant area radius r: 

V^r = -^[(7 + K)at + (7-K)9,*] . (145) 
4k7 

The rescaled vectorfields 

R = , ^ Vr and G = -^^=\/^r (146) 

satisfy the orthonormality relations 

5(i?,i?) = l and g{G,G) = -l and g{R,G) = 0. (147) 

Let g be the affine parameter along R and r the affine parameter along G. In 
the following, we frequently refer to figure [3] of section [3l At the point A we have 
t = T, = VM + ^/2tad by definition. We would like to estimate the value of t at 
the point D and compare it to 1, which is the value of t if f, = 0, T, = \/M and 
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the coordinates are defined on initial data. The rate at which t changes in affine 
parameter along the integral curve of Vj^r going through A is given by 

ar 4k7 — ^ 

We will integrate (|148|) from r = to r = tad with initial condition Tq = \/M + 
y/i-^^ "'^•^ Proposition 14.11 and 14.71 the estimates 

|K + 7-l|<C(e) and , ^ - V2 < C (e) (149) 

hold along the curve. Given the fixed fo we choose the initial data so small that 
To • C (e) is as small as we may wish. Hence tad ■ C (e) is small for any tad (t,) 
with f, < foE3 With these choices the estimate 

\t{D)-VM\<C{e) (150) 

simply follows from integrating (|148p . 

In a completely analogous fashion, by considering and using that I7 — k| < 
C (e) we can show that the point r* = on the initial data is close to r = 2\fM: 
\r{Uata,^)-2^fM\<C(e)^^M■ 

Before we finally estimate how t changes along the integral curve of R through 
D (i.e. the location of the initial data), we derive a rough estimate for the relation 
of r and r* . Consider the vectorfield 

L = ^dr^ , (151) 

whose integral curves are the curves of constant t. The coordinate r* changes along 
such a curve (affine parameter I) by 

dr* _ 1 _ 1 



Integrating from r* = 0, where r w 2v AI outwards to infinity noting that 1 — /i > g 
and that both k and 7 are close to ^ in the region under consideration, we obtain 

q q 

^l<r*<^l. (153) 

10 - - 4V2 

On the other hand, the area radius changes according to 



leading to the estimate 

2VM -C{e)y/M +-1 <r <2^/M + C{e)^/M (155) 
5 10 

Combining (|153p and (|155p yields the relation 

—r~ci<r*<^^r with ci ^ — ( 2Vm - C (e) ^/m) (156) 
11 ~ ~ 8^/2 11 V ^ ' J 



^*Note that f, is close to tad, since the curves = AMp and = iviA converge to one 
ajiother for the initial data going to zero 
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along any curve of constant time in the region r > 2 v AI. In particular, if a quantity 
decays in r in the asymptotic region, it decays in r* as well. 

Finally, we can consider the integral curve of R through D on which the initial 
data is defined. We want to prove that the value of t does not change much along 
that curve (at least up to the area radius R where the support ends). First we 
show that the horizon is a finite length of affine parameter along Vr away from D. 
Namely, since the r-componcnt of the vectorfield R is given by i?'' = \/l — /i, we 
have the equation 

-^v/W. (157) 

Starting at r (0) = l^rtiA and integrating inwards to the point where the curve 
intersects the horizon we find 




< (40HI+C(<5))( 4s-C(<5)) (158) 



^/2 

for some small S. On the other hand, we can integrate outwards from D along Vr 
to a point where r = R. From p55p we know that the affine parameter is controlled 
by the r value along the curve, hence for large R 

g<-R. (159) 
5 

Finally, t changes along the curve according to 

^ = T^(7--)-^<0. (160) 

Within [rK T^y/mA] and [2^ymA, R] we can use the pointwise bound 

{j-n)^^<C{e) (161) 

following from the results of section |4] and choose e (hence the initial data) so small 
that C (e) exceeds the support radius R: 

\tD~t\<C{e)^R<Cie) . (162) 
5 

In this way we can make the difference in t small in the region between r = tk and 
r = i? on the Vr integral curve. 

The pointwise bound of statement 3 follows directly from Proposition 14.11 to- 
gether with the fact that the quantity v is finite in the region under consideration. 

For (|126p of statement V we observe that on {t = Tq} fl {r* > r^} we have 
dtv = by definition. From dr*r ~ (k + 7) (1 — ji) we derive, using Propositions 
I4.1l and [4.7[ the estimate 1 of statement V out = Tq for an arbitrary good constant 
by a suitable smallness assumption on the data. However, along a curve of constant 
r* > r^, the value of r changes only by an amount which can be made small by 
suitable choice of initial data, as is seen from the estimate 



\ritb,r*)-r{ta,r*)\^ / {\ + v)dt< (1 - /i) (ac - 7) < C (e) • Tq 

ta J2VA7 



To 



(163) 
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is small for any ta, tb e 2VM, Tq if the data is small enough. 

The second bootstrap assumption has been dealt with in statement 2 of Propo- 
sition 17.21 already. 

The third bootstrap assumption involves integrals over compact intervals with 
the integrand containing B and its derivatives. The integral is small on t = 
by assumption ^ and Cauchy stability. Again from Cauchy stability it follows that 
will stay as small as we may wish up to the chosen T = d (fo) slice if we only 
chose the data small enough. This is perhaps most easily seen directly from the fact 
that u and v are always finite in the region of integration, and taking into account 
the pointwise bounds on ^,9 established in Proposition 14.11 Put together it 
follows that the quantity can be made smaller than 5 for a finite t slice by an 
appropriate assumption on the data. 

The bootstrap assumptions involving the energy can be satisfied by choosing 
the data sufficiently small (recall the a- priori bound on the mass fluctuation (|45|) '). 
Finally, assumption (|132p follows from the pointwise bound on (cf. Proposition 
14. ip and realizing that integrating the quantity v in u yields a finite result. Hence, 
in the coordinate system defined by tq, all inequalities in the statement V can be 
brought to hold with constant 5 in the region ^. □ 

Corollary 7.1. The set A defined in \133\) is non-empty. 



Proof. By statement 2 of Proposition 17.21 for any f, < tq the coordinates of a 
point in the associated region A {i!) (t,)) will be close to the coordinates of the same 
point in the coordinate system defined by tq. Hence the statement V holds with 
constant 6 in A{d{f,)) for all f, < fg by choosing 6 small enough. Therefore 

[Vm, ^o] cA. □ 

In order to be useful in conjunction with the bootstrap, statement 3 of Proposi- 
tion [7]2] has to hold in any coordinate system Cf associated to a f > fg with f G A. 
The argument is postponed to Proposition 18. 8[ after we have derived appropriate 
decay bounds from the bootstrap assumptions in the next section. 

Proposition 17.21 also provides us with two sources of smallness. In particular it 
justifies the following algebra for constants: 

C{rti)v = 6 (164) 

- ~S (165) 
to ^ ' 

Namely, after we have chosen ?/' > (cf. (|135p and (|136p ) to determine r*j, we can 
choose rj so small that it "beats" any constant depending on r*i, and finally to so 

large that -^^^^j^^ is as small as we may wish. (Of course, the restrictions on the 
initial data get stronger and stronger in this process.) Consequently, everywhere 
that the formulation "we choose t^ so large that" is used in the paper, we always 
have an application of Proposition l7.2l in mind. 



8 Analyzing the bootstrap assumptions 

In this section we are going to derive certain decay bounds for the energy, the 
squashing field and some other quantities. These estimates will be useful for late 
times, i.e. they are to be understood in conjunction with Proposition 17.21 where we 
can choose such a late time. The time to up to which Cauchy stability holds is chosen 
in particular so large that for t > to we have v ^ t in the region r*j < r* < j^t 

and that v ^ t ^ r* in the region r* > j^t. Moreover vo — to + r*i » \fM. 
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All statements about decay in this section are then valid in the subregion {t > 
to} n {u > vq} of the bootstrap region. 

8.1 Energy decay 

From assumption (|129p . we can directly derive ^ decay of the energy in certain 
regions for late times. 



Proposition 8.1. On a hypersurface of constant t we have the bounds 

^ / (-2.) 1!— -^sV^dr* < (t) , (166) 



2t^2 pt-uo 



M 

and 

2^2 pt-uo 



{~2i^) —B^r^dr* < E'^ {t) (167) 



M 



(^{u + af (duBf + {v-af {d^Bfy^dr* <2E§ (t) . (168) 



Proof. The first two bounds follow directly from (|124p . For the last inequality note 
that 2 {duBf {u + af + 2 {dyBf {v - af = (SBf + {SBf and 



(SBY + {SBY = (1 + 2j/) i^iSBY + {SBY ) + (-2^.) ( {SBy + (SB) 
< (1 + 2iy) ({SBf + {SB f 



+4,-2.)((sb + HI;-^s)%(sB + |b " 
+3 (-2,/) (irl^^B' + ^bA (169) 

and that we control all the terms on the right hand side separately by p24p . □ 

The following proposition is an immediate application of the latter and allows 
us to estimate the energy flux through certain slices for late times. 



Proposition 8.2. Let {r\,ti), {f\,ti) be such that ti—f^+a > V M andti+r*—a> 
a/aT and let additionally r* > r"^. Then we have 



m{fl,ti)-m{r*,ti)<3M\ {ti - r*, + a) ^ E'j^ (ti) + {ti + r*, - a) ^ E^ {ti) 
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Proof. 



TO (f^, ti) — m (r*, ii) = J dr*mdr* = J — dum + dyirn^dr* 



7 \ 



where we have used ((30)) . (jSTj) and Proposition l8.ll as well as the bounds ([63]) and 

mi. □ 



The previous proposition can be combined with the bootstrap assumptions (jl3ip 
and l|130p . The fact that energy is conserved then immediately yields decay for any 
achronal slice in a certain subregion of A (T) as elaborated in the following 

Proposition 8.3. In the bootstrap-region A (T) the energy flux through any achronal 
surface 

^C^(T)n{r* < ^0 (170) 
with V- = min„ S > to + r*i > \fM and mint S > to satisfies 

E{S)<M^^^. (171) 



t — t,- 1 1 = i.it,- 



r = 



r — —t 



Figure 6: Energy decay from K. 
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Proof. Dyadically decompose the region A (T) into regions with t^+i = 

l.ltjlffl Proposition IS . 21 applied to any slice ti+i with = r*j^ and f\ = jyti+i yields 
(for late times, i.e. when t — r*i + a ^ t, which is the case for t > to, cf. Proposition 

m (t,+i,r* = ^U+ij -m(i,+i,r^) < Af'^^ (172) 



Combining this decay in the central region with the energy decay at the horizon 
(bootstrap assumptions (|130[) and (|131[) ) we find from energy conservation that the 
energy must decay like through any achronal slice in the region where r* < 

j^t. This shows ()17ip . noting that for large times t > to we have t ^ v in the 

10 



region r*; < r* < j^t. Note in particular that we have this decay of energy flux 



'[ti,t,+i] 

is depicted). □ 



through the regions ^'^"'V^^''^ V '^^ for large ti (cf. Figure [51 where such a region 



8.2 Decay estimates for k and 7 

The following proposition establishes appropriate decay bounds on k and 7 sufficient 
to improve the estimate (jl26p for the relation between r* and r in the central region 
in the next section. 

Proposition 8.4. In the region A{T) n {?'* < ?'*;} n {u > vq} we have 



M 

<Cl^. (173) 



In the region A (T) n {r* > r^.} n {r* < ^t} we have 

\<^<\ + Cl{2 + c)^ (174) 

\>l>\^C,c'-i (175) 
with Ck = sup^*>^*^ j-z^r 0'''^^ Cl ~ sup^*>^*^ 

Proof Integrating equation ([34|) from the set {u = T - r* (T, rx)} U (^{t = T} n 

{r*(T, rx) < f* < r*i}j, where k = i by definition, to any point in the region 
< r*i yields after inserting bootstrap assumption (|132p 

«(i,r:,)-i| <2Ci^ (176) 

in that region establishing ()173p . We can obtain k at any point in the remaining 
region A{T) n {r* > r*j} n {r* < j^t} by integrating from the set L = |{t = 

T} n {r* > U {r* ~ r*i} on which either k is equal to | or satisfies the 

estimate (|176p . to the desired point. An application of Proposition 18.31 then yields 

This decomposition implies that the width of each region is of the size of the t coordinate it 
is at. It should be noted that this decomposition may not fit exactly, i.e. the last of these dyadic 
tubes may have a smaller width. To keep the notation reasonably clean this fact is always to 
be understood implicitly. The results derived for each dyadic region in the paper are of course 
independent of the fact that the last region may be smaller. 
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P7i)l in the region A (T) n {r^, < r* < ^t} as follows 



K{t,r*) 



= K(uL,f)exp — / — — (u,v)du 

< K (ui, w) exp ^sup 
M 



< 



(177) 



For the estimate (|175p . we first note that 7 = i on {t = T} O {r > r^}- On the 
r* = curve we can obtain p75)l by integrating from = T} n {^T < r* < 
T — uq} downwards to any point in the region A{T) n {r* > -^t}. We use that A 

is bounded below and \9\ < the integration region, both following from 

Proposition [131 to obtain 

C (e) M 



7> 



(178) 



there. Since r is controlled by r* (cf. equation ([156])) and r* > -^t in the region 
under consideration, we find the bound (|175[) in the region A{T) n {r* > ^t}, 
in particular on the r* = r^t-curve. Finally, the value of 7 at any point in the 
remaining region A (T) n {r^ < r* < -^t} can be obtained by integrating ([55]) in v 

from some point of the set L' = |{t = T} n {r* < jgilj U {r* — j^t} (on which 

7 already satisfies (|175p ). Using the decay of the energy flux we arrive at (|175p in 
the remaining region: 



7(i,r*) = 7('u,WL')exp - 



> 7(7/, UL')cxp 



It+r 

sup 



r2 A 
2 



(m, v) dv 



r2 (1 - ^) 



[u, v) dv 



t+r* 



> 



1 - 2Ck c- 



M 



(179) 
□ 



From the proof of the 7-estimate we deduce: 



Corollary 8.1. In the region r* > r*j, the estimate ^175^ holds with the constant 
Ck replaced by Cl- 

In the asymptotic region t is like r and the bounds extend; 

Corollary 8.2. In At H {r* > f^t} D {v > vq} we have 



1<^<1 + 2Cl{2 + c)^ 



and 



1 1 

- > T > - 

2 - ' - 2 



Cie) 



(180) 



Proof. The bound for 7 is the statement of (|178p . To obtain the bound for k 
integrate (pi)) from r* = to the asymptotic region of A (T) in u using that 
t ^ r* ^ r in the region r* > -^t (cf. again (|156p ) and that the energy estimate 
holds in the region under consideration. Note again that r could be replaced by t 
in that region. □ 
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8.3 Stability of the coordinate systems 



8.3.1 The relation between r* and r. 

We arc now in a position to derive an estimate for the relation between the co- 
ordinate r* = i^-^ and the function r{u,v). This estimate in conjunction with 
Proposition 17.21 will automatically improve bootstrap assumption [1] which - mod- 
ulo the error-term - expresses precisely the relation of the tortoise coordinate r* 
to the area radius in the five-dimensional Schwarzschild metric. For this section 
we will use C{{rfc, c) to denote a constant which depends on the weight of on 
r = tk and on the parameter c in the bootstrap assumptions. 



Proposition 8.5. The estimate 



Ma (. (r{t,r*)-^mU 
) + \ — [log' 



2 V V^(*,^*) + V2M 



<C[rK,c)^ (181) 



with p defined in {121^ and Ma the Hawking mass at the point (T, r* = 0), holds in 
the region A (T) n {r* > r^}. 

Proof. The estimates 

\dtr\^\\ + v\<C{rK,c)^, (182) 

dr^r = X~^u<{l~ii) + C{rK,c)^< (^1 - + C (r^, c) ^ (183) 

in the region A {T) n {r* > r^} H {r* < -^t} are a direct consequence of Proposition 
18.41 Since also r (T, r* = 0) = 2y/MA, the relation (|18ip follows in the region 
A (T) n {r* > r^} n {r* < ^t}. An application of Corollary [H21 finally extends the 
bound to the remaining region, r* > j^t. □ 

This means that in the region A{T) n {r* > r^} we can go back and forth 
from r to r* with an error-term of j, which is small at late times. In analogy with 
Corollary [HH] we also have 

Corollary 8.3. In the region A{T) {r* > r*j} the estimate U81]) holds with 
constant C (rd, c) replacing C {rK, c). 



8.3.2 Stability of constant t slices 

In this section we are going to study the relation of the different coordinate systems 
Cf = {uf,Vf) associated with different f G A (cf. Section [3]). Instead of the small- 
ness estimates entering the proof of Proposition 17.21 we will now exploit the decay 
estimates for the quantities k and 7 derived in Proposition [831 Recall Notation 13. II 

Proposition 8.6. Let ta G A. In view of Corollary \ 7.1\ assume ta > tq- Then 
in the coordinate system Cf^ = (uf^,Wf^) associated to ta the t-coordinate of the 
initial data slice satisfies the bound 

sup \t - Vm| < C (e) . (184) 

Sn{r*>r*^}n{u<ua} 

Proof. By Proposition 17. 21 the statement (|184p already holds up to tq by a suitable 
smallness assumption on the initial data in all coordinate systems Cf with f < fo. 
Consider now a coordinate system Cf^ for a. ta > tq- Recall the vectorfields G and 
R defined in (|146p . In the following, we again frequently refer to figure [3] of section 
[3l At the point A we have t'^-^ = T by definition. We would like to estimate the 
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value ^Q-* at the point D and compare it to vM, which is the value of i if f = 
and the coordinates are defined on the initial data. The rate at which t changes 
in afhne parameter r along the integral curve of V ±r is given by p48p . We first 
integrate (|148p along the curve = Am a, from A to the point A' , which is defined 
to be on the Vr slice associated with tq. Using the decay estimates 



\k + j-1\<Cl{2 + c)^ 



and 



f 



1 



<C{c)^ (185) 



which hold along the curve by Proposition 18.41 (and its Corollaries), we obtain an 
estimate 

T-&< ^ (ta -rA-) + C (c, ) ^ (186) 



1 - 



To 



where the last term is small and the constant C{c,r*i) depends on the weight of 
on r*j. Using the definition of T = (ta) = VM H — 7== we derive 



4 - 

To 



TA' 



1 



To 



and with the bootstrap assumption on the energy 

TA- 



1 



M 

''fo 



(187) 



(188) 



In the second step we integrate (|148|) from A' to D. In this region we can use the 
smallness estimates for k and 7 as in the proof of Proposition 17. 21 obtaining 



To 



fTA 
''0 



< 



(rA') + VMC(e) 



(189) 



1 - 



where we used the fact -^^C (e) • ta' is small by a suitable choice of the initial data, 

which in turn follows from the smallness of (e) tq by Proposition 17.21 and the 

estimate |t^' — tqI < \ImC (e). Putting together the estimates ()188p and p89p we 
obtain 



'M- 



ta' 



fTA 

''0 



< 



TA' 



'■fo 



MC (e) (190) 



from which it follows (choosing fo large enough and the initial data suitably small) 
that the t-coordinatc at D is close to VM. In the second step, which is identical to 
the one in Proposition [721 one finally shows that t only changes by C (e) along the 
Vr-curve through D on which the initial data is defined. □ 

Corollary 8.4. Bootstrap assumption 2 is improved. 

One easily generalizes the previous proposition to the statement that t does not 
change much along a Vr integral curve located in the bootstrap region: 



Proposition 8.7. With the assumptions of Provosition \8Tm. the t-coordinate along 
the Vr integral curve associated to VM < fi <ta satisfies 



sup \t 

(Vr)^^n{r*>rJ^} 



'd{h)\<C{e) 



(191) 



in the coordinate system Cf^ — {uf^ , ) . 
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Proof. Repeat the proof of the previous Proposition, now integrating only up to the 
Vr integral curve associated with fi. In the second step, when integrating equation 
(jl60|) along the Vr integral curve, one again uses the smallness estimate for k — 7 
in [r'^,R*]. However the decay estimate 

1 M3 

h-n)^=<~e— (192) 



following from Proposition (|8.2p can now be used in the region [R, 00). The e arises 
because is small in _R, 00^ . Inserting that the affine parameter g is proportional 
to r (cf. ([155])), one concludes that 

dt Aft 

< — - < e— (193) 
dg gi 

and hence the change in t along any Vr integral curve is also small within the region 
[7?, 00). □ 



Proposition 8.8. Statements of Proposition \ 77^ holds in any coordinate system 
Cf for f > tq and f A. 

Proof. Bootstrap assumption 1 and the previous Proposition implies that the lo- 
cation of the region "'"'P'^^l^ only changes slightly between the different co- 

[2v M,To\ 

ordinate systems. In particular, the v coordinate of the region ^^^yjj j is 
uniformly bounded in the different coordinate systems, as is the t coordinate for 
r* > r^. Hence if statement 3 of Proposition 17.21 holds in the coordinate system 
Cfo it also holds in the coordinate system Cf for f > f and f G □ 



Finally, we conclude from Proposition 18.51 
Corollary 8.5. Bootstrap assumption 1 is improved. 

Proof. We apply Proposition 17.21 i.e. we choose t^ large such that C (rK,c) 7^ is 
very small (in particular smaller than |) and the initial data so small that the 
bootstrap assumptions hold with constant | at t = tg. Then for t > to the estimate 
of Proposition 18. 51 takes over and improves the constant c in (|126p . □ 



8.4 Pointwise bounds 

In this subsection we derive pointwise decay bounds on the squashing field B and 
its derivatives, as well as on some higher order quantities. The key idea is that these 
bounds hold up to some large time to by Cauchy stability (cf. Propositions 17.21 and 
[8H)E!1 After that time the energy decay derived from the bootstrap assumptions 
in Proposition 18.31 ensures appropriate decay estimates for the fields. 

8.4.1 The squashing field and its derivatives 
Proposition 8.9. The pointwise bound 

\B{t,r*)\<^C{c)^ (194) 

holds everywhere in A (T) H {r*; < r* < j^t}. 

^^The <5 of Proposition 17.21 may have to be chosen slightly smaller but the change is uniform in 

and hence the size of the bootstrap region! 
^''The location oi t = to might change slightly from coordinate system to coordinate system but 
the change is uniformly controlled by C (e) as has just been established in section [8.3.21 
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Proof. The estimate (|194p holds for t G [2-\/M, to] (for some large but finite to) by 
Proposition l8.8l with an appropriate choice of the initial data. For [to, T] we integrate 
out in the u-direction from the set L = {u = uq} U [{t = to} n {r* > jQto)}, where 
either i? = by the assumption of compact support or the bound (|194p holds by 
Cauchy stability, to the r* = j^t curve: 



19 



10 



B[uL,v^—t]+ B^u [u,v ^ —t] du (195) 
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10 



< S 



TO* 4A 



10* -Aki/ 



du<C (e) 



since r ~ r* ~ t on the curve. Note also that along a line of constant v in the 
region r* > -^t we have v ^ t. Integrating out further from a point (t, on the 
r* ~ T^t-curve along the slice t ~ const we obtain 



\B\ < C(e)^+/ \dr*B\dr* 



{dr*Bfr^dr\ 



d 1 

dr* r"^ 



2#r 

or* , 



dr* 



< C (e) — ^ + C (c) —j- sup , 



) dr 



(196) 



yields (|194p in the whole region rj; < r* < j^t since sup . / j,,. < |Cl in the 



region r 



□ 



Recall that in the region r* > j^t we were able to derive -^-decay of the field B 
without involving the bootstrap assumptions (cf. Corollarv l4.3p . The next Propo- 
sition shows that the boundedness of the quantity improves this decay consid- 
erably: 



Proposition 8.10. In the region X ^ A{T)r\ {r* > 3VA/} H {u > VM} we have 

M 



\B\<Cic) 



3 1 



(197) 



Proof. Choose a point (i, r*) in the region X and a point (t, f*) in the central region 
{r*i <f* < 3/M). We have 



r^B^ (t^ r*) = r^B'^ (<, f*) + j Or* (r^B^) dr" . 

J r* 

By Proposition 18. 91 \B{t,r*) \ < Moreover, by Cauchy-Schwarz 



(198) 



^ dr* {r^B'^) dr* <2^ {dr* Bf r^u'^dr*^^ t^B^r^-^dr* 

+ - / t^B^r^dr* . (199) 
t- J ft 

We can finally insert the inequalities (|168p and (|167p to find 



5' (i,r*) < Cl C(c) 



A/2 AP 
2C(c)^ + C(c)-^. 



(200) 
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Noting that in the region it > 1 we have for large times t > ^ yields the desired 
result. □ 



For the region region Uq < u < \fM we can follow the same proof replacing v? 
by + M (to avoid dividing by zero) to obtain 



Proposition 8.11. In the region X = A (T) n {r* > 3vM} we have 



A/ 3 

\B\<C{c) — . 

f 2 



(201) 



Having established better decay of B from the bootstrap assumptions we can 
also derive better decay of 6 via an auxiliary quantity 0, which is the analogue of 
the almost Riemann invariant in four dimensions. Note however that we cannot use 
Q to improve the decay in B itself (as in the four-dimensional case) but only in its 
derivatives, once better decay in B has already been established. 



Lemma 8.1. On {r* = jKt}, the quantity = + ^^\B satisfi 



es 



\Qiu,v)\<C{c) 



■>-'+ 



Proof. Integrate the equation (recall definition ((36|) ) 



35Xv 



2rt 



1 ip2 (B) 
iVF B 



2^/f^ 



(202) 



(203) 



from the set L = {u ^ ua} U {{t = to} H {?■* > -^to)}, where either 9 = by the 
assumption of compact support or the bound (|202p holds by Proposition 17.21 with 
constant 6, to the r* = j^t curve. Since the right hand side of equation (|203p 
satisfies 



B 



35Ai 



4:Vr 2ri 



1 02 (^2 (B) / 3 



4.A^ B 



(204) 



in the region r* > ^t following in turn from the decay of B derived in Proposition 
18.101 we obtain the estimate 



~M3 Mi Mi 

\e {u, v) I < 6 +C{c) < C (c) . 

V-\. r v+ 



(205) 
□ 



Corollary 8.6. 



holds in r* > j^t. 



Mi 

'{u,v)\<C{c) — 



Proof. Use Lemma I5TT] and take into account Proposition lS.lOl 
Proposition 8.12. The pointwise bound 



(206) 



□ 



Aft 

'{t,r*)\<C{c) — 



(207) 



holds everywhere in A (T) H {r*; < r* < j^t} 
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Proof. By the previous corollary 



<(u,v)\<C{c)^ (208) 



holds on r* = j^t. We can integrate equation ([7^ from that curve or t = t, 



0, 



where the bound \6\ < holds by Proposition 17.21 to any point in the region 



A{T)n{r*<^t}: 



{u,v) = e{u,,v)~^ / ^Xdu+ / _ (e-8B _ e-2B) (209) 



2 Jui f Jui 

and hence 



. X , 3 / r 4C2A , / r -Aki^X , 

'(.,.)! < ^—^-Ji^d^\ i^d^ (210) 



The energy estimate, Proposition 18.31 and the fact that v ^ t in the region r*^ < 



< Yr;t yields the desired result. □ 



With the previous Proposition and Proposition 14.11 Corollary 18.61 is easily ex- 
tended to the entire bootstrap region: 

Corollary 8.7. 

\0{u,v)\<Cic)^ (211) 

holds in all of A (T). 

Close to the horizon we have 
Proposition 8.13. The pointwise hounds 

\B {t, r*i) I + \e {t, rl,) I < 2 C (c) ^ . (212) 

hold everywhere in A (T) n {r* < r*,} n {w < T - r* (T, r^)}. 

Proof. The decay for 9 was already obtained in the proof of Proposition l8.12l From 
Proposition 18. 91 we know that on r* = r*i wc have the decay 

\B\<^C{c)^ (213) 



Consequently, 



B{u,v) = B{ur*^,v)+ / —{u,v)du 



< vci:c{c)^+ . r ^du. r -^du (214) 



and upon inserting bootstrap assumption (|132p we obtain the result. □ 
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8.4.2 Higher order quantities 

In this subsection various bounds on tlie derivatives of tlie quantity fl^ are proven. 
Since we iiave not yet establislied a pointwise bound on ^, tlie estimates will turn 
out to be suboptimallf^ However, they sulRce to estimate certain error-terms in 
the X-vectorficld identity. An interplay between the X and the y-vectorfield will 
finally generate a pointwise bound on ^ , which allows one to optimize the estimates 
(cf. Proposition 18. 16|) . In particular, the new decay will then suffice to control the 
error-terms occurring in the identity (jlOp for the vectorficld K . 

The first step is to improve Proposition 14.61 to a decay bound. In the following 
C {r*i , c) denotes a constant whose weight is determined by Cl and which also 
depends on the c in the bootstrap assumptions. 

It should be emphasized again that the quantity is only piecewise contin- 
uous, with a discontinuity spreading along the null line v ^ T + r* {T^tk)- The 
estimates below are valid because the quantity du-^, which is integrated along 
null-lines, is continuous. The same considerations are valid for the quantity 
whose discontinuity is along the null line u = T — r* (T, rx ) . 

Proposition 8.14. In the region A{T) n {r* > rj^-} n {r* < ^t} we have the 
one-sided bound 

^-;:,<C(.:.,c)^. (215) 
In (T) n {r* < r*,} we have 

-^--,<C{r^,,c)-^ (216) 

Proof. Since k = i on the u ~ T — r* (T, r/^) ray we have 



/4k 


4 


K 














" 3 


r 





Q^2K-^-K\^m+^'-\p-"-]] (217) 



and, in view of Proposition [8T3l the estimate 



^{u^T-r*{T,rK),v) < !^ (u ^ T - r* {T,rK) ,v) + C (r^^, c) ^ 



Moreover, on {t = T} n {r* (T, rx) < r* < -^T} we have by the constancy of k 

= 2.^-4^-+-- p-- -.-^ (219) 



and hence 



fly m , ^ , \fM 

-^--,<C{rl,A^ (220) 



following from the fact that \B\ < — in that region by Proposition 18.131 Note 
that the inequality (|220p would also be two-sided if we had the analogous pointwise 
bound on ^. Integrating ([26]) downwards from the set L = {u ^ T — r* (T, r^)} U 
{{t = T} n{r* (T, r^) < r* < r* }) to the r* = r* curve yields 



— iu,v) = — iu,,v) 



K.mv ^ V ( 3\ C^' 
-6 — ; 2k— p — - — 3' 



j,2 y 2 y 7-3 



{u,v)du (221) 



^^Such a pointwise bound could in principle be established via a bootstrap argument in the style 
of Proposition 14.21 with the pointwise decay bound on B l|213| l now entering the estimates. 
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and upon inserting the pointwise estimates on B and 9 fProposition l8.l5|) , bootstrap 
assumption [5] and the estimate 



-3 



CO 



(u, v) du < 3C (c) 



f72 



du < C (c) \/Cl — 5 — 2" 



,2 V4 



for which (|132p has been used, we finaUy find that 



^ <C(rJ;,c)^ 



(222) 



holds everywhere in r* < r*j estabhshing (|216p . Starting from this curve or from 
the curve {t = T} (1 {r*^ < r* < j^t} we can integrate further to any point 
in the region A (T) n {r > r^} H {r* < j^t}, this time using the energy estimate 
instead of to obtain (PT5)) . 

□ 

Proposition 8.15. In the region A (T) n {r* > r^} n {r* < ^t} 

Cie) 



< 



(223) 



the region A (T) n {r* < r^.} n {u < T - r^} 



> ^ > - ^''^^ 
Q, ~ r^ v+ 



(224) 



Proof. 7 i on the set L = {t = T} n {r* (T, rn) < ^T}. From 



= 27- 



n 



47 47 



3 r 



-7 --T"^- 



1?1 

V2 A 



(225) 



we derive using the decay estimates (|194p . (|212p . (|207p the bound 



on L. We write the evolution equation ([^5]) as 



A'/ 



(226) 



A 2A 



9. ( ^ ) - 7 ( 6m-j + _ ( p- - ) + 3---3 



C A 



(227) 



and integrate downwards in v. Using the estimates (|174p . (|175p and again the decay 
estimates for B and 9, the error-terms are estimated 





37 T 




J V 







M4 

(u, v) dv < C (c) sup 



and 



7— 



< 



C(e) 



(228) 
(229) 



This establishes the estimate (|223p in a subregion (w > -^T) of the region asserted 
in the Proposition. For the remaining part, we derive the estimate 



771 



(230) 
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valid on {t = T}r\{r* > ^T} using the decay of B, 9 { CoroUarv IS . 71 and Proposition 
IS.lip in r. Integrating (|227p downwards to any point in the region {r* > j^t} using 
again the estimates for B and 6 one obtains (|230p in the entire region {r* > -^t}. 



Since t ~ r* on the curve r* 



9_ 
10 

n 



t we obtain 



<C{e) 



M 
rt2 



(231) 



on that curve. FinaUy, integrating (|227p from the r* = j^t curve downwards up to 
any point in the region r* > yields (|223p for the entire region asserted in the 
Proposition. 

For the estimate (|224p we integrate (|227p from r* = where t ^ v downwards. 
Clearly the round bracket on the right hand side of (|227p is always positive and 
hence the upper bound of ()224p follows immediately. For the lower bound we use 
the estimate 7 < ^ available in the region under consideration to estimate: 



A 2A 



- 6m— + -t[p-- 



1 



1 



> ^ (u, M + 2r^) - m (m, u + 2r^) {u, v) ^ (m, u + 2r^) - 



lit 

> ^ - 



We easily extend the bounds to the asymptotic region: 
Corollary 8.8. In the region A{T) n {r* > j^t} we have 



□ 



(232) 



Proof. The first bound follows from integrating equation (j26p outwards from r* = 
jlyt using that r ~ t in that region and the decay of the fields in r. The second 
bound was obtained in ([230ll . □ 

As seen in the proof of Propositions 18.141 and 18.151 a pointwise decay bound 
on the quantity ^ would considerably improve the estimates on the higher order 
quantities. We summarize this as 

Proposition 8.16. Assume that 



< C holds in r* < r*i 



< C 



Ml 



in —t>r*> r*, 



(233) 

holds for a constant C depending only on r*j. Then in the region {r* > r*} for any 
'i^ci ^ ''g ^ ''if '"'6 have the bounds 



m 



<C{r*i,c) 



t^ 



and 



m 



M 



t2 



(234) 



Moreover, the one-sided bound 121 6\) in the region A (T) n {r* < r*;} is extended to 



IT 



<C{r:i,c) 



(235) 
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and the bound 



is refined to 



0>—{u,v)> — - {u, v)-C {r*i,c) 
the region A (T) n {r* < r* } r\{u<T - r^}. 



(236) 



Proof. Revisit the proof of Propositions 18.141 and 18.151 Note that the constant in 

the M-estimate of (|234p improves by moving away from the horizon since it depends 

on the weight of on r* □ 
° 1— u g 



9 The vector field Y 

RecaU the functions a and f3 defined in section [71 Close to the horizon we are going 
to apply the vector field 



for which 



2a 



Y = ^^^du + 2P{r*)d, 



= 2/3 y„ = -(3n^ n = -a . 



(237) 



(238) 



The calculations will be carried out in the Eddington Finkelstein coordinates defined 
in section [21 From ([5^ we derive the identity 



1 { 2 



2 {daB) 



I , av ^ l3Xn'^ , 1^, 
2 



f]2 



il'^r \4k / 

In a characteristic rectangle TZ = [ui,U2] x bi,i'2] the identity 

F^{{U2} X [vi,V2])+F]; {[ui,U2] X {V2}) 

= ll (TZ) + Fl {{ui} X [vi,V2]) + Fj ([«i, U2] X {vi}) 
follows, with the boundary-terms given by 

±^Fl{{u}yi[v,,V2])^2 ' (^P{d.Bf + ^[l-'^f^c^r^dv, (242) 
1 



Fl{[u,,U2]x{v})^2 



and the spacetime-tcrm 
1 



V2 rU2 



V2 rU2 



Vi J Ui 



J_/_2 \/l f3Xn^ _ 1 ^, 

'27-2 I 3'^/ r r 2 



(239) 



(240) 
(241) 



^ + (l - ^ ) ) ^''^^ (243) 



^ ( 4a— - a ) + [dyB) 



(244) 
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It will be useful to split the term into 

/^(7^) = -/^(7^) + /^(7^) (245) 

where 



1 /•V2 /•«; 



f72 V VL 



r^dudv (246) 



ana 

—/J = y _/ (- [j^'^'^^) duBd^By^dudv. (247) 

With the choices of the functions a and /3 made in section [71 the integral (7?.) is 
non-negative for r* < Vy and moreover, using (jl38p . 

il (?■* <r^) = 2tt^ J^' J ' (^-^a - duBdySy^dudv 



V2 



^ £" ^ (^"^ - + ^' 

< 5-fl(r'<'-y) (248) 

holds in r* < r^. We conclude by rewriting identity ()240p for a characteristic 
rectangle with one boundary being the horizon: 

{{Uhoz} X [vuV2])+F]^ {[ui,Uhoz] X {w2})+/b (7^) 

= /^(7^)+FJ({^il} X [vi,V2]) + F^ {[ui,Uhoz] x {vj) . (249) 

10 The vectorfield X 

All calculations in this section are performed in the Eddington Finkelstein coordi- 
nate system defined in section |3l 

10.1 The basic identity 

The vector field X is defined as 

X ^2f{r'')du~2f{r^)d, (250) 

for some function / chosen below and with uj satisfying uj > ti — r*j. It will be 
applied in the region 

■="H=t2-r:,jyr:„uj .251) 

[11,121 [tl,t2j ^ ' 

for some r*^ also chosen below. We note 

X" = 2/ X" = -2/ Xu = fn"^ Xy = -fVt^ . (252) 
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From now on primes will denote a derivative with respect to r* , hence d^f (r*) = 5/' 
and duf {r*) = -i/'. From ^ using 



(253) 



we compute 



x-^ {n^) -(n^) \ 

-r-^f+ ^ ^'" / .(254) 



With the boundary terms 
1 



+ 



{duB f (2/) 



fdtBdr*B{U,r*)r^dr* 
2 



and 



^ ttX 

27r2 



27r2"^" 



X 



vi 

2t2~Uj - 



l--p]{-2f) 



du 



/ 2ti ~uj 

one can state the identity 



{d^BY (-2/) + _ ( 1 - -p ) (2/) 
(9.5)^-2/) + ^ (2/) 



(255) 

(256) 
(257) 



/ ^ [-T^.^^-- ~ (V^T^,) X^] dVol = (ti) - F# (to) + - Jl . (258) 



[tl.*2l 



Let us turn to the spacetimc integral on the left of (|258|) with the integrand being 
given by (|254p . In view of definition p6p we can write 



W-BW^B^laB' + ^{e' 
and the integrand (|254p becomes 



5-'^) = ^nB'-8^ + ^¥'2(i?) (259) 



T^uTt'^-' ~ (V^Tps) X' = —r {dr^Bf + -UB' (_/'__ (A - ^) / 



1 



+ -^(^1 (i?)) -/ 



7' 

^2 [B) 



-/ + 



-/'--(A-i.)/ . (260) 



Finally, we apply Green's theorem to the □i?2-termP^ Collecting the i32-terms of 
P60|l after the integration by parts we find 



-52 



+ n{f' + -{x-u)f 



(261) 



^^cf. the remarks in appendix 1X1 
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Since 



Dw{u,v) = ( --^dudy - ^-^dy - ^-^du ]w{u,v) 



(262) 



and moreover / depends only on ?'* we arrive at 



f' + f 



6 / X-iy\ 9(A-i^)" 



/ 



12 n o 



18 ^ 



A- 



18 A „ 



A-z/ 



Computing the derivatives explicitly for the expression in the square brackets of 
(p6T|) yields the identity 



- 32 



12 



dudv 



18 V 

To? 



18 A 



du 



^2 



35 18m\ J_ a. 



(35 + 9^) +I7 (5,0,0 (263) 



with 

It{B) 



_9_^ ^ 36 AC^ 



16 



3\ 1 



2/ r3 



-14(A- I/) +8r 



-2B „-8B 



)(^?-C) 



(264) 



We summarize the remaining error-terms as 



- V ( ipi {B) + 3(^2 (-B) A \ , (/^i {B) + (^2 (S) 



f72 



52 



52^2 



and read off the pointwise estimate (cf. Corollarv l4.3p 



X7 {B) + Xg (S) < C (e) 



M 



(265) 



Taking care of the boundary terms arising from the application of Green's identity 
we can finally state the identity (cf. equation piip ) 



X jX 



(266) 



with the renormalized bulk-term 



52 
2 



r 6 fx-u 

172 r I ^2 



r + r 



6 

?72 



[*l>t2l 

\-iy 



9 {X-vY 



r2 f]2 



y r2 V f7 

1 
2 



B'- 



1*1, *2l 



l7iB)+Is{B) 



>dVol 



dVol. 



(267) 



53 



the new boundary terms 



' [r + -{X-i^)f] (dtB) B {t, r*) r^dr^dAs. 



S3 



- at - (A - i.) / W?^ (ii,r*) r'dr^dAs^ 



/' + ^ (A - / ) (duB) B{u,t + r*) r^dudAss 



+ 



t-r* JS 



r 



the horizon terms 



J^UH - ^UH 



Bd^B [J' + -{X-iy)f 



and the J-terms 



- .F 



du[ f + ^{X-iy)] f] B^{u,t + r*i) r^dudAs. , 

(268) 

{uhoz,v) dv 
r^Uho..v)dv (269) 

(uj, v) dv 
r^{uj,v)dv. (270) 



2t2—UJ 



2ti-u 
2*2 -ii J 



Bd^B [f + ^{X-u)f 



2ti —uj 



-ir^v /' + - (A - z.) / 



10.2 Analysing the X-bulk-term 
10.2.1 Borrowing from the derivative-term 

We would hke the spacetime term (|267p to have a sign. To achieve this we borrow 
from the term containing a derivative. Define 



, h + 12 



* „* , *1 ~ *2 



and rt = r 

Ki +ti-t for ti<t<t' 
r*i+t- h for t' <t<t2. 



(271) 
(272) 
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and computi 



40 



dt 



t—u.i fl 



f72 



dt 



t-Uj rl 



{dr^B + £,Bfr^n^dr* 



dt 



dt 



Kit) 



) ^ ^ ~^ '-^ error ^^ Horror 



X,U2=t2 — r*, 



for some function ^ chosen in (|277p . The boundary-terms are 



fiB'r' {t,t-uj)dt 



X,vi=ti+r* 



dtf^B'^r^ {t, VI - t) 



and 



H. 



X,U2=t2 — r*i 



dtf'^B^r^t,t-U2) 



2t2-Mj 

2ti—uj 



f'^B'r'iuj,v)dv, (273) 



f'^B^r^u,vi)du (274) 



f'S^B^r^ {U2,v)dv . (275) 



To keep the notation clean we write M = m{T, r* = 0) in this section. For a 
sufficiently large constant a we define the shifted coordinate x 



We choose 



X ^ r* — a — V M . 
3 X — v nx 



2 r x^ 



(276) 
(277) 



for some n £ (ijCX)) from which 



e' + e'--(A-^)e 



9{X-iyf x^ {n^ - n) + na'^ 



4 r2 (a2+a;2)2 

follows. Hence the integral (|267p can be expressed as 

2f 



\^ltl,t2] J 



ri al-"-' 02 

'tl.*2l 



{dr^B + ^By dVol 



- + f 
02 



4n.T 



02 ((j2 + a;2) 



+ / 



(ct2+x2)2 02 



+/ 



35 18aA 1 /O^^, O 



(^--) -i:^--f +-r Tf-^ (35 + 9M) 



o o 



B2 dl/o/ 



("^^7(s)+X8(s)^s^dyo/ 

tA',«j , „X,ti=ti+r;, „X,u=t2-r* 
error ' terror i 



(278) 



"'"Again care is necessary in the integration by parts because of the difFcrcntiability of the 
coordinate system. In any case it is sufficient to note that the integrands of the boundary terms 
are continuous and that the integrand of the bulk term is picccwise continuous (all terms except 
the ^'-term are in fact continuous everywhere). 
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which wc will write shorthand as 



(-p,rli,uj\ _ jX ('r)'-li.uj\ jX,uj 



where 



TjX,v—ti-{-r*i jjX.u—t2 — r* 
Tiprrnr ~r terror 



F + f-g--IjiB)+IsiB) 



(279) 



^dVo{280) 



with the identifications 



1 



F = — — — f" + + f 



35 18// , I f 



(281) 
(282) 



10.2.2 The choice of / 

By Propositions 18.41 18.141 and 18.151 we have the bounds 

M 



n n - r t 



< 



1 



'M 



C{e) 



in the region r^j < r* < j^t. This implies that 



and that 



9- lyf (35r^ - 104Mr2 - 10871/^) 



< 



C(e) 
tr^ 



Mr- 



(283) 

(284) 
(285) 

(286) 
(287) 



where M = m{T,r* ~0). It is apparent (note Proposition 14.61 in particular) that 
the expression g is negative close to the horizon, positive far away from it with a 
single zero at some r = Vzero (t) for which the estimate 



2A/y ^ (26 + Vl62lj < C (e) Vm 



(288) 



can be derived. This in turn implies an estimate (with error of order j) for the r* 
value along that curve via identification p81|) . It follows that g changes sign in a 
small interval around 7'J„^„. For future calculation the estimate 



(289) 



for will be sufficient. 



We finally construct the function f = f (x) = f [r* — a ~ V M \ by setting 



a/ - r* 



(290) 
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and 



M"-2 

(.t2 + (72)" 



(291) 



Note that / will be bounded for n e (^00). Later we will set n = |. We compute 



-2xn 



f" = AP 



1 2ri (x2 + 2nx-2 - cr^N 



(a;2 + 0-2) 



to find 



, - 0-2 



f]2 (^2+^2)"+2 



(292) 
(293) 

(294) 



We will now show that there exists a positive constant c (a) > such that in 

F + f-g~\lAB)+I^{B)>^c{a) (295) 
2 r-' 

holds in the region of integration. Note that the above choice of / ensures that 
f ■ g is positive everywhere, except for an e-small correction-term. In particular, in 



r* -\/lVf r* + — \//V? 

' zero 10 V , ' zej-o ^ 10 ^ 



we have 



fg>-—TC{e) (296) 
and outside of this set, using (|265p and (|286p 

l/g-ix7(S)+l8(B)>^. (297) 
10.2.3 Estimating /f (^^li^'t^^ ) 

The aim is to establish (|295p . We will do the computations in the shifted x- 
coordinatc (|276p . 



The region x < —a and x > a Note that F is already non-negative for x < —a 



and X > a. In the subinterval 



where / • g might cause problems, we have the stronger bound 



F>-^ci (a)>0, 
M 2 



which upon combination with 



yields 



the only subset 
(298) 

(299) 



for the regions under consideration. 
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The region [— cr, a] Wc shall show that wc can dominate the term F in the region 
[— CT, cr] by the term |/ • g in (|280|1 f3 In conjunction with (|297|) this will yield the 
desired result (|295p . By Proposition [ 



- (1 - /i) I = I (47K - 1) (1 - /i) I < C (r;,, c) ^ (300) 

holds in the region — <t < x < cr and also / is positive therc|3 In view of (|286p it 
suffices to show 

M"--^^^^ < A f35 - 104^ - 108^) /I (301) 
(1 - m) (a2 + 2;2)'+" - 2r3 V r2 r4 ^ 8 ^ ' 

in the region — cr < x < a. There wc can estimate 

/ (X) = r f > r /' > M-i . (302) 

J-ff-\/I7+r;„„ J~a~^ {2a^ + 2<tVM + M] 

such that it suffices to establish 



(303) 

Part /. Consider first the region — cr < a; < — |cr translating to \/M < r* < 
icr + ^/M. The lower bound on r* implies a lower bound on r by identification 
()181|) . In particular, r > \/M (2 + |) in that region and hence 

1-M>T^ and -^<^ (304) 
10 1 — /i 7 

as well as 

35-^-108^ > 17 (305) 

hold in the region under consideration. Consequently, for cr sufficiently large we 
have 

n a"^ — 10 cr — X a + x 

< n- 



(1 - fi) (^2 + x2)2+" 7 V^^T^ (fa2)'/'+" 



-10 / 9 ^ 



The upper bound on r* can be exploited for large cr to give 



r<^<7. (307) 



Again choosing cr sufficiently large, this gives rise to the estimate 
1 A35_ 104M_^pg_A/2\ x + a + VM 7 



2cr2 + 2crVM + M 



17 /12y ,r+^7 116 x + a 7 
- 2 ^5<Ty (3a2)"8- 3«a2»+3 8' ^ ' 



^'^Note that we save g/g to get the positivity of II297I I. 

*^Note that we can use the aforementioned Proposition because the region under consideration 

lies in r* > r*. . 

— cZ 
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Comparing (|308p and p06p for n = | we have shown the desired inequahty (|303p 
in the region under consideration. 

Part II. Consider now the region — fcr < x < a, which translates to + \fM < 
r* < 2a + \/aI. We can choose a so large that 

1 - M > y (309) 



and 



' <(l) (310) 



X + a \6 

hold in the region under consideration. We deduce that for large a 



7 1 /^^ 104m ,„„TO^\ x + a + VM ^ 7^2 x + a 6 
35 7^ 108— r rrr > :r35- 



8 2r3 
and 



{2a^ + 2aVM + My ~ « 5 (2^^)" (x + a)^ 7 



< . (311) 

— / o . '->\2+n ^ 



(l-Ai) (a2+x2)2+" - 6 (^2 + ^2) 
Hence it suffices to show that, for large a, we have the inequality 



3 



1 „(a-xHx + a) 

9 ^ ^ (a2+x2)'+" 



For n I we obtain 

lV2(--)(- + f -\l, (313) 

3 (CT2+a;2)2 

which is shown to be true by elementary arguments. Namely, for x < we have 
l_^^{o^x)ix^of^ < 1^2^^^ < V2V2F^^a < 2 ^ , ^^^^^ 

3 (cr2+2;2)2 3 (^2^^2)2 3 (^2_|_2,2)2 3 

and for a; > we have 

V2 (cr - a;) (a; + af ^ \/2 cr* (x + af ^ ^2 ( ax + 



3 (ct^+x^)^ 3 (0-2^ ,^,2-) 2 3 Va;2+a-2/ 

This finally establishes that the integrand of is non-negative and vanishes if and 
only if _B = everywhere. 

Remark: The minimum size of the constant a required for the estimates above 
to work can be determined explicitly. Choosing a = A0\/M for instance is large 
enough. 

10.2.4 Summary 

We can write 

I ^[ti.ta] ) ~ {'-^[tiM] ) ^error J^error — terror l-J-^Oj 

and we have shown 
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Proposition 10.1. The estimate 

l,...^dVol<Cia)I-{vlt:r^ (317) 



[tl.'2l 



holds. 



In the following subsection we are going to estimate the error boundary terms 
on the right hand side of (|316p . 

10.3 Controlling the error-boundary-terms 
Lemma 10.1. The error-terms (273^-(27S^ satisfy 

\Jerror\ < C [a] (m (uj, 2^2 - uj) - m [uj, 2tl - Uj)) (318) 

HfrrZ=*'^'*' >0 (319) 

and 

HfrroT'"'"^^' > . (320) 

Proof. The statement pi8p is immediate since we are away from the horizon and 
both ^ and /' decay sufficiently fast at infinity to retrieve the correct powers of r 
appearing in the energy. For the other two inequalities recall that /' > always 
and that x < —\/M in the region of integration and hence ^ > 0. □ 

10.4 Controlling the boundary-terms of (^[fi tz] 

The following Lemmata show that the boundary terms of the vectorfield X appear- 
ing in the vector field identity (|266p are controlled by the energy plus a contribution 
from the term appearing as bootstrap assumption [SI Together with the results 
of Lemma 110.11 identity (|316p ultimately yields an estimate for the positive space- 
time integral manifest in Proposition ll0.2l 

Lemma 10.2. We have, for any q E IR+, 

it) I < [C (r^,) + Cfq] Ep (t) + CfrlB^ {t, r,*,) + -CfPl {v ^ t + r*i) 

(321) 

with 

Ep (t) =m{t,t- uj) - m (uH, t + r*i) (322) 

and Eg (wi) being the quantity appearing as bootstrap assumption Moreover 
Cf = supr<r*^ [f'f] <C 1 is a small constant and C {r*i) just depends on r^^. 

Proof. The i^-boundary-terms arising from energy (|119p can be estimated 



..3^(5.«r 2 



{u,t + rli)du 



^r'-^iB,,r + Ar'—iB.ur + riX-,y) ( 1 - §P ) ) it,r*)dr* 



2 



> 



{u,t + r*i)du 

t-r!, L V - / J 

-J ^ I {dtBf + {dr^Bf + —\r^t,r*)dr*{?>2?>) 



Cl 
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The F-boundary terms arising from the basic identity for the vector field X (|268p 
are 

f^t — Uj 



fdtBdr*B{U,r*)r^dr* 



(324) 



f + -{X-iy)f]{dtB)B{t,r*)r^dr* 



r^duBfi2f)+'^(l-lp) i-2f) 



t2 — r* 



du 



/' + - (A - z^) / ) [duB) B{u,t + r*) r^du 



7;\du[f' + -{\-v)f]]B^ {u,t + r*i) r^du 



Spacelike i^"^-terms We estimate the first three terms of (|324p (with the index 
si denoting the restriction to the spacehke terms) 



(t) I 



t—uj jr 



{dtBfVM^ ^ 



M , 

t — U J 



r^dr* + 



< 



* I/I {{dtBf + {dr^Bf) r'dr* 



t — U J 



\f\^{X^„)[r{dtBr + —]r'dr 



I/I 



2r^ 



{X-v)i\ + iy) + - (r^„„ - r „„) 



B^r^dr* (325) 



Recah that /' is positive and that we have |/| < C (a) and f'<C (a) Moreover 

(A + i^) is clearly bounded everywhere, as is and r^%^ (cf. CoroUarv 18. 8p . 

Using the equations 



2^X-le^ and r. 



O 9 

2 a2 



,C (326) 
it becomes apparent that the terms in (j325p are controlled by energy. 



Null F'^-terms For the last three terms of ([324]), which only arise for the (ti)- 
term by the geometry of the region, we observe that the first is manifestly controlled 
by the energy because / is bounded along the rays of integration. For the second 
term one estimates for any q G 



f\duB\\B\THu < Cf 

ti—r* 



qrB^-u) + -^-^^r^ 



du 



where bootstrap assumption (|132p has been used, and 

'3 



(327) 



< 



/3 



Iti-r* 



-(A-^)|./|U 



1 / [duB 



(A-i.)|/| \duB\\B\{u,t + r*)r^du 



M + -^B^ {~v) (u, vi) du<C (a) Ep • 
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Finally, for the third term we note that x < — VXT and hence /" > in the 
integration region to estimate 



2 



-duf - {duf) - (A - z.) + I/I du[-{X-v) 



2|B||B ^Ir^dw + C (m (wi, - m {u2,vi)) 



t-[ —r*. 



1 ^, 



< Cfrii B' {h,r2i) + Cfq ■ Ep + Cf-Fj^ (i-i) + CEp 



for any q G 

Lemma 10.3. We have 



\Jl\<CEj{uj) 



for some constant C and 

Ej (uj) = m {uj, 2t2 - uj) - m {u,j, 2ti - uj) 



□ 

(328) 

(329) 



Proof. For the terms in the first line of (|270p apply the inequality 2BB y < + 
r (B^v)^ to retrieve the correct powers of r in the energy. For the second line observe 
that /', r_yy and r_uu decay sufBciently fast in r. □ 

Lemma 10.4. We have 

<C Eh {vi,V2) + CfrlB^ [t^, K^) 

for some constant C and 

Eh ivi,V2) = m {uhoz,V2) - m {uhoz,vi) . 



(330) 



(331) 

Proof. The term is manifestly controlled by the energy by the fact that / is 
bounded. For the second and third term in (|269p . we observe that the terms which 
are multiplied by a A- or a v-i&cior (or derivatives thereof) are controlled by the 
energy. For the remaining terms we note that x < —\/M and hence /" > in this 
region and estimate (using a Hardy inequality) 



B^ 



\B\\d..B\f' + — {f)^^ 



r^dv 



<GfrtiB' {t2,rli) + 



B^fr^Xdv 



2\B\\d,B\fdv 



< CfrliB^ {t2,rli) + CfEn + 16 H {d.Bf ^-^r^dv + / 



B^ 



(/'),„ 



r^diii2) 



from which we obtain 



B^ 



because 



|B||a„B|/' + — (/'),„ 



dv<Cfrt,B'{t2,rli) + EH (333) 



M 



(334) 

is small in the region under consideration allowing us to estimate the derivative 
term in the last line of (|332[) by the energy. □ 



62 



From the previous Lemmata and the identity (|316p we conclude 
Proposition 10.2. In the region 'D^f^'^J^ we have for any q G M+ 

If (ijg;,"! ) < [C (r^,) + Cfq] [Ef ih) + EF{h)]+C (6) ^ 

+ CEH{vi,V2) + CEj{uj) + 2Cf^(^F^ ivi)y (335) 

Proof. Add up the estimates of the previous Lemmata. Note that pointwise bounds 
for on the r* = r*j-curve were obtained from the energy Ef (t) in Proposition 
18.91 including a small term C (e) Hence we can replace the terms appearing in 
Lemmata 110.21 and 110.41 by the energy, such that (|335[) is finally obtained. □ 

10.5 Controlling the time derivative from (B) 



For the statements of the next two propositions let i?* = ^ 5 V M < T*zero ^'^^ define 
the regions 

B^{h<t<tl}C^ {rli <r* < R*} (336) 
and the slightly smaller region 

= {to < < < <i} n {r*; +^<r* <R*-^}. (337) 

where <r < is some positive number (in units of a/M), say <^ = j^VM. 

Define also a smooth cut off function x supported in [r*i , R*] and equal to 1 for 
{r*i + < r* < R* — Note that with the definition of r*^ (cf. section we 
have r*i + (; < ry and also R* — > — i-y/M. 

Proposition 10.3. We have, for large ti, 

^ B'dVol < C (a) I§ (Sg;;;) (338) 

and 

^ (dr* Bf dVol < C {r^i , a) /f (sg f [ ) . (339) 
Furthermore, we can control the time derivative 

-±=l^id,BfdVol < C(r*,a,x)[/f 

+ m{t2,R*)-m{t2,r*i) + m{ti,R*)-m{ti,r*i){3A0) 

Proof. Inequality (|338p is the statement of Proposition llO.il Equation p39p follows 
from an application of the triangle inequality to the first term in (j280|) and the 
previous bound on the B^-integral, noting that ^ is bounded above and below in 
the region under consideration. Finally, (|340p is obtained via Green's identit^^^ 



.t = t2 

□X ( -is'] = j^xi-\^B^^ - J xBdtBr'dr^dAs. 



(341) 



*^There ajre no problems with the differentiability here. 
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for the X defined above, which can be written 



{dtBf dVol 



V2 (B) 



dVol 



J X {dr^B f ^dVol + j xBdtBr^dr*dAs3 



The spacetime integrals on the riglit hand side are controlled by p38p and p39p . 
For the boundary term in the second line we estimate 



xBdtB{t,r*)f^df*dAs3 < 



B^ 



+ VM [dtBY ] {t, f*) f^dfdA^z 



< VMCir*,)im{t,R*)-r7iit,r*,)) . (342) 
Putting all this together we obtain 

/ {dtBf dVol <^ [ x^ (dtBf dVol 
VM Jr.- vM Jb 



We can summarize this as 
Proposition 10.4. The quantity 

1 



□ 



w 



(dtB) 



Lidr^Bf + ^B^ 
M Mi 



dVol 



(343) 



satisfies 

Ib{B,) < C(rJ„a,x) [If 

+ m (<2, i?*) - m (t2, r*i) + m {ti,R*) - m (ti, r*;) 
Proof. This is a direct consequence of Proposition 110.31 



(344) 
□ 



11 Combining X and Y: Horizon Estimates 

For the next two propositions recall the choice of Vy made in section [5J which 
implied in particular that Y is supported only in r* < — i\/M . 

11.1 Controlling from Jf and energy 

Proposition 11.1. Consider the characteristic rectangle TZ ~ [ui,Uhoz] x [vi,V2\ 
together with the r* = r*i curve intersecting {ui, vi). Define u {v2) by r {u{v2), V2) — 
r*i and r(u,v{u)) = r*j. We have the estimates 

(uuvi)) (345) 
:}) (346) 

) for all u>ui (347) 



Fb X [vi,V2]) < C{r*i){m{ui,V2) -m 

J U(V2) 



Ul,Uhoz 



X 



v{u) 



rB'^dv < CF]^ (^{u} X vi,v{u) 
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Proof. This follows from the definitions ((M^ and □ 

Proposition 11.2. Recall the basic dyadic regions 2?^^"^' j for the vectorfield X 
i251\) and erect the characteristic rectangle 

n = [ui,Uhoz] X [vi,V2] (348) 
associated with such a region as depicted in the figure. More precisely, let ui = 




Figure 7: The horizon estimate. 



Let also 

T = {r* > r^-} n {w < V2} n {it > Ml} , 
and recall that R* ~ —^\/M. We have the inequality 



(349) 



<Cir*i,a) 



{B[tTM ]+r>iit2,R*)~^T^it2,r*i)+m{t,,R*)-miti,r*i) 



- j ~r iiL 1^1,^, J-i. ) — III, yu2, I clJ ~r 

+C{r*i) m{ui,V2) - m{ui,vi) + Fg {[ui,Uhoz] x {vi}) .{350) 
Proof Recall the identity ((249)) : 



Fb {{Uhoz} X [vuV2])+F}^ {[ui,ui,oz] X {v2})+iB {T^) 
= il (7^) + Fl {{Ui] X [VI,V2\) + Fl {[ui,Uhoz] X {vi}) . 

By Proposition 1 11 . 1 1 we control 

Fb ({"1} X [vi,V2]) < C {r*i) (m (ui, U2) - m (ui, wi)) 
To establish pSOp we will show 

1 



(351) 
(352) 



il{n)<-il{n\T) 



+C {r*i,a) [/f (Sg'j) + m {t2,R*) - m {t2,r*i) + m {t,,R*) - m {t,,r*,) ] (353) 
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+m {t2, R*) - m (t2, Ki) + m {h^R*) - m (ti, r* ) (354) 
To see this decompose 

il{H) = il{n\T) + il{T) (355) 

Since \iiTZ\T we have by definition r* < one can apply (|248p to obtain 

/I(7^\r)<i/^(7^\r). (356) 

On the other hand, in the region T we have 

lliT) < C{r*y)lB(Tn{r* <R* = -^^}j , (357) 

which follows from the fact that Y is only supported for ?-* < — i\/]\7. An appli- 
cation of Proposition ll0.4l to the term on the right hand side of ()357|) will produce 
the required second term on the right hand side of (|353p . The estimate p54p is 
obtained completely analog ous to ((357)) . □ 

Proposition 11.3. With assumptions and geometry as in Proposition \11.2\ we also 
have 

F^{{uhoz} X [vi,V2]) + F^ {[ui,Uhoz] X {w2}) + i/^(7^\r) 
< C {r*i,<j) m ^ = ti - R*, V = ^t2 + R*^ - m{u = t2- r^i^v = ti + r*;) 

Proof. Use the inequality 

/f (sg;J)</-(— (358) 

which is obvious from the positivity of the integrand (compare the dashed lines in 
the previous figure for the regions). Inserting the estimate of Proposition ll0.2l into 
the inequality p50|) we obtain the result by an appropriate choice of q. □ 

It is of crucial importance that the constant C {r*i) just depends on the choice 
of r*i and not on rj^-. 

11.2 Controlling from and energy, on a good slice 

Finally, we are going to control the boundary terms by and the energy on a 
"good" null-slice. 

Proposition 11.4. With TZ and T as before pick a v e [wi,U2] that satisfies 

{[ui,Uhoz] X {v}) = inf Fj {[ui,Uhoz] x {v}) . (359) 

Vl<V<V2 

Then 

Fl ([mi, Uhoz] X {v}) <C{V2- viY^ il [n \ T) + C ir*i) (m (ui, W2) - m (iti, t-i)) . 

(360) 
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Proof. Recall that the expression (|243p is manifestly positive. Set u{v) = w — 2r*j 
and estimate 



{[ui,Uhoz] ^ {v}) < inf Fg {[u{v) ,Uhoz] ^ {v}) + Fb {[ui,u{v)] X {v}) 

VI <V<V2 



< 



V2 ~ Vi 



Fg {[u (v) ,Uhoz] X {v})dv + C {r*i) [m (wi,f2) - m (wi,fi)] 



where v is the ti-slice determined by taking the infimum of Fg in the region 
[u (v) , Uhoz]- For the integrand of the first term in the last line we have 



F^ {[v - 2r:i,Uhoz] X {v}) < 2n^ / du rHVM 

'v~2r 



r r 2 it 



,(361) 



following from the fact that the inequalities (|139p . (jl40p . (|14ip hold in r* < r*i. 
Comparing (|36ip with (|246p produces the first term in (|360p . □ 

We will also need a related version of the previous Proposition, which provides 
one with a good energy-slice instead of a good F^-slice: 



Proposition 11.5. With TZ and T as before pick a ij £ [vi,V2\ that satisfies 



E {[ui,Uhoz] 'X {v}) ^ inf E {[ui,Uhoz] x {v}) 



(362) 



Then 



E{[ui,u,,oz]x {v}) <C{v2~vi) ^il{n\r) + {m{ui,V2)-m{ui,vi)) (363) 



Proof. Recall that 

E{[ui,Uhoz] X {v}) 



^ (9„i?) 



is manifestly positive. With u (v) ~ v ~ 2r*j estimate 
E{[ui,Uhoz] X {v}) < inf £;([u('(;),Mfto2] X {«}) + £;([«!, u(u)] X {w}) 

Ul <V<V2 



< 



E {[u (v) , Uhoz] X {v}) dv + [m (wi, 1)2) - m (wi, wi)] 



The integrand of the first term in the last line can be controlled by 

^2 



E{[v-2r*i,Uhoz] X W) < 2tt 



2 / M, ^3 



du r-^WM 



(Mr f 4.^ - a' 



2r' 



(364) 



following from the fact that the inequalities (|139p . ()140p . (|14ip hold in r*j. Com- 
paring (|364p with (|246p produces the first term in p63p . □ 

The results of this section are already sufficient to obtain a pointwise decay 
bound for the quantity ^. For reasons of presentation this is postponed to section 
113.11 but the reader impatient to see the argument can turn to the latter section at 
this point. 
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12 The vectorfield K 



12.1 The basic identity 

The vectorfield K is defined as 



(365) 



It is the analogue of the Morawetz vector field in four dimensions. In particular, it 
is conformally Killing in five dimensional Minkowski space|f^ We note 



M 



9? 2 

{v - a) Kv 



M 



(u + a) 
"(366) 



and 



u ~ t — r* V ^ t + r* (?j — a)^ — {u + a)^ = At (r* — a) . (367) 



From ([5^ we compute the identity 

M (-T^.^^'' - y''%sK^) = ^[i^iu + af + X {v - a)') OB^ 



+32^ - ^ (" + «)' + M« - af) + ^ ((1^^) „ (u + af + {n') {v - af) 
+ ^ (" + a) V A (v - af) (ifii (B) + ip2 (B)) - ^ipi (B) (i/ {u + af + X {v - a)^) 



(368) 



with ifi and (y52 defined in (|36p . We shall apply the basic vectorfield identity in the 
region (ef. Figure [5]) 



for any T < T producing the identity 

n (V-^^,,) = F- (f) - F- it,) + H 



^ ^ . +0 

UH=1 —rl. 



where 



1 

M 



1 

M 



[to,Tl 



t— Mo 



(369) 
(370) 
(371) 



[to>Tl 



{duBf 2{u + af + {d„Bf 2(v- af 



(372) 



+ ({u + af + {v-af)^{l-\p 



r^dr* 



2 (li + a)' r3 (a„B)' + _ (^, _ a)^ ( i _ 



and 



27r2 Af 



T+r* 



2 (t; - a)' [d.Bf + [u + aYil- '-p 



(?i, t + r}()du 



T — r , 7J I (if . 



^''The vectorfield field has been shifted by a for reasons which will become apparent later. 
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Note that the J-tcrni vanishes in view of the assumption of compact support. We 
are now going to define the renormahzed quantities Ig and Eg that arise from an 
apphcation of Green's theorem to the OB'^ term in the spacetimc integral (|37ip . 
The D in the basic identity pOSp is here given by (cf. appendix 



D 



3 f ly {u + of' + A (w — of' 



(373) 



We compute 



2 



V llL + of' + A (w — of' 



= t -24r- 



A A 
12r— r,„„ + Urr^uv— - 24— (r,„„),„ + 24 



+ + a)^ [ [A + H(^ + |^) -§((-,„),. + (....),..) 



+ {v- ay 



3A 



3i/ 



(374) 



and define the bulk term 
1 /■ /■ 1 



+t 



2 



32 - 24r^-12— + 4^1^ 



A^i' JTutjA r"^ , , A 
64A + 24— + 12^^ - 24— (r „„) „ + 12r— r,„„ 



f72 
^1 (i?) 



-8r- 



17. 



(u + aY 



\ n 

QvX 16 1 ^1 {B) + 3v?2 (^) 



4A + 12A 



^2 (i3) ' 

B2 



^2 



+ i'T.uv)^ 



16 + 2- 



n ^ n ) ^172^,™ 3^2^,' 



(375) 



In order for the identity (cf. equation (|llip ) 

to hold, the boundary terms have to be 
(t) 1 rt-^" 



T-rt 



(376) 



27r2 



27r2 



+ /. ^"'2B{dtB)(^[,.iu + af + \iv-af))r'it,rndr* 



t-rt 



+ m/ . 2B(a„B) f-(t.(« + a)VA(i;-a)2)jr3(,i,t + rJ,)d7. 



^ ''-B^a^ (^|- (i.(ti + a)" + A(«-a)") ) r^(?/,t + r^)d?/ (377) 
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and 



27r2 



1 ^.fif..2.., .2^, 3 



2B{dyB) [ — {iy{u + ay + X{v-ay] ] r^T - r*j^,v] dv .{378) 



12.2 The Spacetime integral 



Let us turn to an analysis of the integral (|375p . Besides the formulae (|326p . 
and (|25p the following identities will be useful 



(?',Ut)) ,j 

Q2 



O „ 2 / 3\ A / 3\ 4 



-2B 



Q u 2 f 3 \ 



r2"3r V 2 



3r2 2 



3r 2 



-2B 



The bulk integral can be written 



{'^\to.f]) ~ ^B,main {'^^„^f]) +^B, error (^^^.T]) 



(379) 



with 



35 + 9/X + 4 



52 



8 p- 



+t 



r ~ a 



24^r^-64r^ + (l-M) 



-70k — 36K/i — 



X380) 



3 



12/i 



+ 



-35 - 18m - 14 ( p - - I + 2 ^5 



+ iv- ay 
where 



16+-(1-m) + 2— ^ 



6^7^ + 6- 



>(381) 



32 



^-2 



-28k (1 -/i) + 16^r 



(382) 
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and 

QiB)='-T- 16 (e-- - e--) ^ . (383) 

Note that 

\P{B)\<Cie)^ (384) 
by the pointwisc bounds of section 

12.2.1 Estimating /f 

We start with the observation that I b main has a good sign near the horizon and 
near infinity: 

Lemma 12.1. One can find R* such that the integrand of I ^ main negative for 
r* > R* . It is also negative for r* < r*j. 

Proof. The second statement is a consequence of p4ip and the inequahty < 
which follows from Proposition l8.l51 For large r* on the other hand, we can expand 
the integrand of (|380p in powers of ^ using the results of section [ 



and (cf. identification (|181[) ') 



!_ ^ (^^j where p = ^ ^ ^ ^ 

to find 

2 



2V2 + log 



2 + 



„2 



With the chosen centre a of the K vector field (a = —p — 1 by equation p27[) ). the 
integrand will be negative in r* > R* for some suitably chosen □ 

Remark: In particular, we will choose to so large that R* < j^to holds. 

The idea in estimating the spacetime integral I § main ij^\t f]} decompose 
the region of integration into dyadic pieces (cf. footnote [ 



^^„.f] = E^fe,*.+.] witht^^f (385) 

^^„t,+,] = 2^^,,T] n {tj < t < tj+i} ■ (386) 

For each piece 2?^ ^ 1] control the bulk term Is^main by the bulk term 

losing a power of t: 



^^Note that the rest terms are all controlled by £iil^-!£l. 
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u 



2v - 4r 



K 



Figure 8: Different regions to control the error-terms of K. 
Proposition 12.1. In the region T)^ we have, for each dyadic piece 

Proof. By the previous Lemma it suffices to show (|387p with 2?^^ replaced by 

' because the integrand of I § main admits a good sign to the left of ?'*; and to 
the right of R*. For the compact r'^-interval the first part of inequality p87p follows 
from Proposition llO.li the second from R* < -^^to and the positivity of □ 

12.2.2 Estimating /f ^rror 

In this subsection we are going to show that the contribution of the integral I b error 
can be made as small as we may wish for late times by choosing r*j^ sufficiently close 
to the horizon and the initial data sufficiently small. To achieve this we shall split 
the integration into different regions U,V, W and Z defined as follows 



U = Vll^^^ n {r* < rj,} n{u>2v~ 4rJ,} , (388) 

V = Pf^^^j n {r* < r*j,} n{u<2v- 4rJ,} , (389) 

W = Vllf^ n {rj, <r*< At} , (390) 

^-^£,.f]nK^^^}- (391) 
An immediate observation is 
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Lemma 12.2. In all regions we have 

/t + 7 > 0. 



(392) 



Proof. This is a consequence of the choice of coordinates and the monotonicity of 
K in M and of 7 in ?; manifest in equations p4p and (j35p . □ 



The next lemma estabhshes appropriate bounds to control the error-terms of 
Ib error cquation (|381|) . 



Lemma 12.3. Recall that by Provositions \8.(A \8.12\ and \8.13\ the bound 



\B\+AI-^\9\<C{r*i) 



t 



holds in W and 



\B\+M-\9\<C{r:i) 



M 



holds inU UV. Assume also 



Mi 

<C(r:,)— inW 



<C(r*,) in UUV. 



Then we have the following estimates 
• In region W 



\QiB)\ 



6 



M4 



7* 2 



t2 



<C{r),,c) 



M 



t2 



• /n region V 



6- 



< 



7-2 172 

C(r*,c)^/I7 



> 0. 



• /n region hi 



v < di exp — 



/or positive constants c?i > 0, c?2 > 0. 
• /n region Z 



Proof. 



(393) 
(394) 

(395) 



(396) 
(397) 

(398) 
(399) 

(400) 



(401) 



The region W The bound (|397p is the statement of Proposition l8.16l The bound 
P96p follows directly from the decay properties p93p and p94p . 
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The region V. From Proposition 18. 161 we derive the bound 



by observing that u is like v in the region V. 

The quantity is obtained by integrating (|26p written as 



(402) 



A 2A 



= 7(6?7i— + — (p-- ) +3-- — 



OCX 



(403) 



from the set L = {{t = T}C^ {r^. < r* < r*,}} U {r* = r*J downwards. On L itself 
we have by Proposition 18.161 



(404) 



Since 7 < ^ in V by monotonicity and moreover will always be negative, we 



can derive the bound 



n 



{u, Vr) 



A 2A 



7 1 + ^2 



e C A 



m y/ M 

> — ^{u,VR)-C[r^i,c) — :r- + m[u,VR) 



> 



r3 
m 



(m,u) - C (rj,,c) 



i2 



1 



(m, w) (ii; 
1 



(u, vj^) {u, v) ^ 

(405) 



in V, where in the last step we used that in V the Hawking mass decays like 
Putting the boimds (|405)l and ((402)) together yields ((3991) . The bound p98)) on 
Q (B) follows directly from the pointwise bound p94p . 

The region U Integrating the quantity 1/ = 7 (1 — /i) from the spacelike t = T 
curve downwards to any point in the region U we obtain 



("7 v) = ^ (1 - m) ("t, vt) exp 



- (1 - 1^1) {uT, vt) exp I - / / (u, v) {u, v) dv 



with 



~ 4k 4 k 



/ {u, v) {ut, v) dv 
" f{u,v) {u,v)dvj (406) 

(407) 



In both regions V and U the quantity / is clearly positive, bounded below by some 
^2 > 0. We can estimate, for a point (w, v) in region U 

- iy{u,v) <^{1- fj.){uT, Vt) exp J f {u, v) {u, v) di?j 

< ^ (1 - /^) (uT,WT)exp (-da (wr*, - w)) < i (1 - /i) (uT,WT)exp ^-^-^^u^ (408) 

Here we have used that v < Vr*^ — by definition of the region U. 
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The region Z The estimate is tlie statement of Corollary (|8.8p . □ 

With the necessary bounds in place we can prove the following 
Proposition 12.2. Assume iS95\) holds. Then the error-term I § error satisfies 

N-l 

+ Mi{r)i) + M~6{tQ) (409) 



with 



lim e(rjf) = as well as lim 5 {to) = 0. (410) 



Proof. By Lemma 112.21 we have A + > 0. Hence the term multiplying (A + v) in 
Ib error ^ good (negative) sign in all regions and can be ignored. For the other 
two terms we look at the different regions: 

Region W: Note that u and v are controlled by t in this region. We insert (|397p 
and p96p into the integral I b error- The resulting term, which has to be controlled 
is 

VmC / du dv n^r^— . (411) 
Jw 

We split the region of integration into = Wn{r* > r*;} and = W n {r* < 
r*i}. The region is partitioned into dyadic slices as was the bulk term: 

AT-l 

I§,error (W') = E ^^.e™ + ■ (412) 

We can control each dyadic tube by Ig losing a power of t (arising from a missing 
power of r in (|41ip ') 

I^rror K,,*,,.]) < (r J, , a) ^ + i [/f (W[l^.,^.^^] 

< C{rt,.a)J2t,+, 



J=0 



In the region we can ignore the factors of r. It suffices to estimate < '^^'^f 
from Proposition 18.91 and hence 

dtdr^n^r^^< fdt^^l^ r'dr^^^<Cir^,„c)^. (413) 

Region V: In this region u is like v at late times. We insert the bound (|398p into 

I B, error ^"^^ estimate the resulting term [vq = t^ + r*i, V = f + r*^) 



C dvduVL^B^< / / du(~AKv)<C—. (414) 

M Jv Jvo VM^^Jta-r* ' ' Vo 



For the remaining terms, i.e. those containing a + -factor, we insert the 
one-sided bound (|399p to control the error term 

f dtdr*r'^^u^dtn^) (415) 
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for large times (in partieular to + r'^- > 1) as follows {uh = T — rj^). First define 

f-r' 



-J + r^ and 



r_ = Ti 



' K 



and 



tin 



to + Tk ~ f* for r* > 
~3r* + 4r^ for r* < 



(416) 
(417) 



Then the term (|415p can be estimated using Proposition 18. 131 

rf-rt+r* 



1 

M 



dr* 



,5^ 



t(r*) 



dtr^—u^[dtn^ + cn 2 



< ClC (c) Vm / dr* 



dt dtn^ + cn^- 



/M 



Jt(r* 



< ClC (c) %/m / dr*n^ (T-r*ji + r*, r* ) + C^C (c) M / dv du— 

Jr*_ ^ ^ Jta+r* Jv-2r* ^ 



<ClC{c)VM r T,r^ -r T-r^ + r*,r* +ClC(c)M 



dv—;: 



ta+r*. 



duVL' 



/M 



< CiC (c) M • e (rj,) + ClC (c) Af (rj^ ) 

to 



where in the first step we have used that the round bracket in the first line is 
positive. The constant C (r*;) may have different values in each line. We also used 
that 



4 V 7 K 



and therefore 



7K 
7 + K 



< 4a 



(418) 



(419) 



holds. In summary, smallness for this error-term arises from the smallness of the 
r-difference between any two points in the region r* < rj^. The crucial point is 
that only C (r*j) enters the above estimate, such that the r-difference can "beat" 
the constant. 



Region U: To control the error-terms in region U estimate the curly bracket of 

error somc coustant timcs and by something small (cf. Corollary [ 
The resulting integral can be controlled via (|400p as follows: 



\lB,error W I < TT^ (e) / dv du (^, y) 



< — C(e) 



to+ri 



rT-r*^ 

dv I ducxp 

2i.-4rl- 



u]u^ <MC{e)Ce~^'r 



2VM 



The region Z On the one hand, we have to establish smallness for 



1 f 

. dudvr^^f 

M 



{u + of + {v — of 



< CVm J dudv 



Q{B)+Q 
Q{B) + 6 



r202 2Kr2 
AC2 3(?2 



2^2 2Kr^ 



(420) 
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where we used that r controls v and u in the region under consideration and Propo- 
sition 18.111 From Proposition 14.11 it is apparent that the critical term to control 
is 



1 1 



1 



(421) 



M 

Namely, the remaining terms in the square bracket of (|420p all decay like ^ by 
Proposition 14. II such that direct integration will already lead to a smallness factor. 
The other critical term to control is 



1 /■ 1 , 



M 



((m + af + (i; - af ) 



(422) 



which upon inserting (|40ip and using the fact that r controls u and v in the region 
under consideration reduces to controlling the term 



1 



T-'uo 



dt I dr*B'^r (-v) . 



(423) 



Using that r ~ t in region Z we can estimate (|423p by 



< 



1 1 /■-' 

C{r*,,c)^ dt-. dr*B\''{-v) 



< C {r*, ,c)^ I dt^E§ (t) < C {r*i , c) ^ 



Ml 



where we have used bootstrap assumption (jl29[) . 

12.3 The Boundary Terms 

We write the boundary-terms p77|) as 

(0 ~ ^B,mam (0 ^" ^B.errorarc (^) ^" ^B,errorline (^) 



(424) 
□ 



(425) 



where 



B ,main 



(t) 1 



27r2 



M 



t — Mo 



-12i uBdtB + 6B^iy- 



r^dr* 



M .L. 



t — Uo 



{duBf 2 {u + af + {d.,Bf 2{v- af 



M 



l-T-rl 



t-r* 



2 (u + af {duBf + ^ (« - af (l - 



(M,t-|-r*,)dw, (426) 



PB,errorarcii) 1 ^ 2 f i'^ + ^f , , , (u - fl)^ \ 3 



1 ( +3BdtB^^^^^ (A + ^^) - 3^^B^ iX + ^)] r'dr* 

VI J^*^ \ r r J 

B^ (^^^ (i^ {u + af + X{v- af^^ r^dr* (427) 



27r2 M J^. 2 

pt-ua ( (',, _ 

2^ i-t—uo 2 



M L* 2 
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and 



B ,errorline V / 

2^ 



and 



1 

M 
1 

M 



T-r 



f-r 



2B 



27r2 M 



2 (v - af [dyBY + 11^ (u + a)" ( 1 



{T-r*K,v)dv 



T+r] 
ta+r* 



Note that F|;^„„,,,„^ 



^ J" 2i? (^A (y {u + af + \{v- af)^ (T 

(f) 



r^,z;)di(429) 



= and that the last term of F^^^^-^^ also vanishes for 



t = T. 

12.3.1 Estimating F|;,„,,„ (t) 

We are going to show that the boundary term Fg^^^^^^^ (t) comes with a sign. This 
is obviously the case for the integral in u, so it remains to establish non-negativity 
of the spacelike integrals. Define 

S ^{v^a) a„ + {u + a)du, (430) 

S^{v-a) a„ -{u + a)du. (431) 



Note that 
and 

respectively 



(SBf + {SBf = 2 (it + af {duBf + 2{v~ af {d^Bf (432) 



S^tdt + (r* - a) dr* S = tdr^ + (r* - a) dt (433) 



tdt ^ S - {r* ^ a) dr* tdt^ , J , S--^ -dr*. (434) 

[r* — a) (r* — a) 

We can insert these expressions into the boundary term (|426p and integrate the 
second term by parts using S: 



T~uo 



1 



BtdtB- {-2v) {r* - a) r^dr* 



T-uo 



-2v 



[r* -a)B[ {SB) - [r* - a) dr*B]r^dr* 



= vr^ [r* - af 



'-B[SB) 



T*=f-UQ 




4 


/ (-2zy)r3 




hi. 



-B' 



r* — a {r* — af 



dr* 



(435) 
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where 



TV 3 



C(e) 



(436) 



is very small by Proposition l4.1l Note that the boundary term near infinity vanishes 
in view of the assumption of compact support on the initial data and the domain 
of dependence. The term at r* = is manifestly positive since u < Q in the 
integration region. 

Alternatively, using S_, we can obtain 



T~uo 



.1 



T-ua 



-2v 



BtdtB- {-2v) (r* ~ a) r^dr* 
r 



r* —T—uo 



(r* -a)B \ -SB - -dr^B r^dr* = vrH^B^ 

^ ' ^ {r* ~ a)- {r* ~ a) > 



T-uo 



i-2iy) - (SB) B + —B' 



^-[iy + r-^]+PiiB) 



dr* 



Again the boundary term has a positive sign at r* = rj^. 

If we split the relevant term in (j426|) into two equal pieces and collect terms we 
can write 



o_? B .7nain 



1 



M .L* 



T-uo 



{duBf 2{u + af + {d^Bf 2{v~ af 



^^u + af + {v~af)^(l^-p 



r — a 
r 



I2t- -vBdtB + 65^1^- 



r^dr* 



> 



M 



T-uo 



dr*r^ 



(1 + 2iy) (^{SB f + {SBY 



i-2iy)\ [SB + l^^^^B 



2 r 



23 



+ 6 + 3 



K — V + r- 



^.u 



3 t 



2 +2 



23 



{-2v) \ [SB + --B] +-B^[- + 5 + -i 



which is manifestly positive. Furthermore 



K — V + r- 



B ,7nain 



(r) > i?f (f)+ - I {-u) — [e + {r* - af) r'dr* (437) 



with being the quantity appearing in bootstrap assumption [31 



.errorarc 



T 



12.3.2 Estimating 
Lemma 12.4. Under the assumption Ii395\] we have that 



B .errorarc 



T] > [T] + Mi 



with the e arising from the fact that T is large. 
Proof. We have 



\\ + ,y\<C{r*^,c)^ for rl.<r*<Af 



(438) 



(439) 
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by Proposition 18.41 and 



\\ + iy\<C{r*j„c) 



M 



9 ~ 

for r* > — T 
- 10 



(440) 



by Corollarv l8.2l Recalling the bound (|8.16p we can estimate the expression 

3' 



:— ; 

r 



"37 



(441) 



by 



|r ™ + r.„,„| < ^ (A + z.) + C (rj,, c) ^ < C (rj,, c) ^ 



and similarly 



\<C{r*K,c) 



'M 



J^2 



both in the region r*^^ < r* < j^T. Analogously, we obtain 



+ \r.uu + r. 



uv I 111 2 



c 



9 , 

in the region r* > ^g-^- 



(442) 
(443) 

(444) 



Inserting these estimates into (|427p it becomes clear that we have to establish small- 
ness for the terms 



T-Mj 



o o 3 6^ 3 C 



T, r* dr* 



(445) 



We split the integral into the region rj^ < r* < j^^T and the region r* > jqT. In 
the first region the derivative terms of (|445p arc manifestly controlled by the energy, 
decaying like ~ by Proposition 18.31 The term can be estimated as 



B'^r^dr* < Cf / S^rdr* < Cf 



' ' 10 



T 



(446) 



In both cases smallness is obtained from the fact that to is chosen very large. In the 
region r* > j^t on the other hand, the derivative terms in (|445p can be controlled 
by pulling out the i as a smallness factor and use the energy estimate for the rest. 
For the B^ term wc have to borrow from the good last term of (|437p : 



10 



T-uj 

B^r^dr* < 
If 9T 



T~uo 



B^r^dr* 



(447) 



Hence a tiny contribution from the last term of (|437p will control this term and we 
finally arrive at □ 



12.3.3 Estimating H, 



K 

T-r 



Lemma 12.5. Under the assumption 11395]) we have 



(448) 



where e (r|f ) for r*^ 



00. 
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Proof. The first term of (|429p is elearly positive and can be neglected. For the other 

terms spht the integration X = [vi, V2] = to + r^-, T + r^- into the part which lies 

in U (where we can use the estimate (|400p ') and the part in V (where we are going 
to exploit the fact that the r-differenee is small) . See Figure [S] Following this line 
of thought we estimate the (negative of the) second term of (|429p 



rV,^„ {v - ay + 2r^\ (w - a) T - r^, w 



< 



[/ ' 


I ] 


3^2 / 


Jxnu 


JinV. 


2 \ 



+M2C(r*„c) 



< C7e"5"(u + a)^ 



inu 



B^rX {T-r*f.,v]dv + CM^ 



'inv ^ 



dv 
dv 
dv 
dv 
dv 



+M'C{r*,,c)e{r*j,) < M^g(rJ,) (449) 



where we used that r^uv < 0, that u is like v in region V, and the assumptions p94p . 
For the (negative of the) third term we obtain (C just depends on r*j) 



T+r* 



2B (d^B) 



ta+r* 



< 



inu Jinv 



+ r{B,,f^ i-i.) {u + af ( 



T ~ r'^,v] dv 



B^ 

r 



+ l\r^{^+r{B.^fY{v^af(^l 



< Ce" 2" (u + a? + CMi 



inv 



X{T-r*K,v]dv + CM^ /AT- rt,v] dv 



T — t'^,v] dv 



r-r^,w) dv 



<HrKim 



where we again used that u is like v in region V, and the inequality ()394p . as well as 
the fact that (— i^) < A (cf. Lemma ri2.2p . These estimates together yield (|448p . □ 



12.3.4 Estimating {t) 

Clearly ^^^^^^^^^ (^T^ ~ since there is no upper null-boundary for the region in 

which we apply K. Hence we only have to estimate Fb errorHne i^o)- This is done 
in the same manner as for the horizon term: Splitting the integral into a part lying 
in V and a part in U, using the estimate pOSp in the former and applying (|400p in 
the latter region. 



12.4 Summary 

We have shown the following 

Proposition 12.3. Assume ^3951) holds. It follows that 

N-l 
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and e can be made arbitrarily small by choosing both —r*^ and then sufficiently 
large. 

Proof. Write ((575)) as 

B,main [ J B ^errorarc \ j B^raain I [^q T\ j B .error \ [ti.t2] J 

(to)-i^t=f-.;, (451) 

and apply the estimates of Lemmata ll2.4l and ll2.51 as well as Propositions 112.1] and 



13 Closing the bootstrap 

With the required estimates now in place we are in a position to prove the closed- 
part of Theorem 1 7. II i.e. to improve the remaining bootstrap assumptions 
We start with the observation that the X-bulk-term decays. 

Proposition 13.1. We have 

{VI::C"') < Ib {t^:;^:^-') < . (452) 

Proof. Apply Proposition 110.21 in combination with Proposition 18.31 and the boot- 
strap assumption p32p . □ 

With the help of the Propositions proven in section [11] we can derive the point- 
wise bound (|395p , which was assumed for most of the Propositions established in 
section [nPl 



13.1 A pointwise estimate for ^ using Y 

Proposition 13.2. In the region A (T) n {r* < D {v > to + r*i} we have 



V 

9 



MS 

<C(r:,,c)— (453) 



and in A (T) n {r* > r*,} n {r* < rf^t} the estimate 



Ml 

<C{r*,,c) — . (454) 



Proof. Starting from the slice Sj^ (cf. definition 1 116p erect the characteristic rect- 
angle to any 'St, to < t < T. By Cauchy stability (Proposition 17. 2[) . we have that 

i[ui,Uhoz] X {vo =to- r*i}) < S (455) 
and hence an application of Proposition [TTT5] together with (|452p immediately yields 
^Fj {[u,,UHoz] x{v = t- r:,}) <^S + e (456) 



*® Recall that the first two have been improved already in Corollaries 18. 41 and 18.51 

*''The reader is assured that none of the results of section 1121 will be used in the following 

subsection. The argument has been placed in this section because it is also used to improve the 

integral bound I I132I I. 
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for any t <T. This estimate and Proposition II 1 .21 immediately imply the uniform 
estimate 

1 , 11- 

(457) 



for the region TZi\Ti of any characteristic dyadic rectangle. Next we apply Propo- 
sition 111.41 to each dyadic rectangle to find a slice v satisfying 



1 

M 



FU[^^,^ho^] X {v})<C- 



O IT 



(458) 



Proposition 111.21 applied to the rectangle enclosed by the good slice in [vi, w^+i] and 
V = Vi-f-i yields, again using (|452p 



(459) 



Having exported the better decay to all late slices in this fashion, we can erect 
the characteristic rectangle again and apply Proposition 111.21 which produces the 
uniform decay estimate 

1 . /Tr 

(460) 



^/I(7^A7^)<C^ 

M V,, 



One may repeat the procedure, i.e. apply Proposition lll.4l again. which now provides 
one with a good slice (with the F-flux decaying like jj^-^-yz)- After application of 
Proposition 111.21 this leads to the decay 
1 



—Fb {[u;,Uhoz] X {n,+i}) <C{rli,c) 



M 



(461) 



M - — iv,+,r 

on any late slice Vi. Finally one may export the decay to any w-slice by choosing 
appropriate regions: 



1 

M 



F^i[u{r*i),Uho,]x{v})<C{r*i,c) 



M 



1 



M 



— F^(^.x [«,«]) <C(r*,c) -2 

M 



-il{n\T)<C{r%c) . 2 



(462) 

(463) 
(464) 



everywhere. Note that we have assumed a better bound than (|462p in the bootstrap 
assumption (|132p . however the bound (|463p is new and essential to derive the point- 
wise bound for ^. Namely, integrating (|57p upwards in a characteristic rectangle 
yields 

— (u, Wi+i) = — (w, w,;) e ^ iT\i' 2;j^ ■ > 



and hence 



361 4 k 



2r 3 



iu,v)dv (465) 



3 1 



2r. 



min v J Vi 



Vi+l 



c c 

-{u,Vi+i) < -{u,Vi) 
V V 

-/j[ilfH(«.«)'^V(«,^;)dz;^£'"^ 



{u, v) dv 



+C sup 

\ r<r*, 



arB"^ (u, v) dv 



Mi 
<C + 

Vi 



.l-d-Vi 
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Reiterating from the first to any chosen late rectangle we find for any {u,Vi) in the 
region r* < r*^ 



C M3 C Mi 

- iu, v^) < C {r*,,c) and hence - (u, v) < C (r*,,c) (466) 

v Vi y V 

which is ((^551) . Integrating dHZD from the set L = {r* = r* }U{{t = to}n{r* > r^,}}, 
where the bound (|454p holds by Cauchy stability and the estimate just established, 
we obtain (|454p in the complete region using the energy estimate and the fact that 
u ^ t in the region where r*i < r* < -^t. □ 

13.2 Improving assumption (11291) 

With the pointwise bound on ^ established we can improve assumption ()129p for 
any late boundary term Ej^ (^T^ via Proposition 1 12. 31 One applies the JiT-estimate 

in the region "^'"""^^^^^j^' j'^j some large —r^, late to and tjy — T. Using (|452p we 
have 



fX I u=t 



Mi 



(467) 



with the e arising from the fact that to can be chosen as large as we may wish (at 
the cost of making the data smaller). Consequently 



N-l 



N-l 



^ t,+i/f H^+^"'--2?fct'^'"0 =^ ^(*o)A^^ E < e(io)Mi (468) 



[tj .tj+i] 

in view of the finitencss of the geometric series 



'i.r 



(469) 



Tl = 



Combining (|468p with the fact that Fg (to) is small by Cauchy stability, Proposition 
yields 



for any late tN, which improves assumption (jl29p . 



(470) 



13.3 Improving assumptions (fT50]) . (fTMl) and (fT^ 

We apply Proposition 18.31 again, inserting the better bound for the A'-boimdary 
terms (|470p at late times to improve the decay of the energy on any arc-part of late 
shces, St n {ift > r* > r*,}. From Proposition [8Jl we also obtain improved decay 
for the field B in the region r* > r*i. 

Proposition 111.51 applied in each characteristic dyadic rectangle produces after in- 
serting the better energy decay in the region r* > r*i a slice Vi = Ti + r*i with 
improved energy-flux decay, ttt^ + ttti- By the domain of dependence property 
the decay of energy flux is improved on the horizon piece v € [vi,Vi^i] and on the 
ceiling-part of the characteristic region, as indicated in Figure [9] This retrieves in 
particular assumptions (|130p and (|13ip with a better constant. In view of the energy 
flux decaying now like on all achronal slices in the region A (T) fl {r* < jgi}, we 
can apply Proposition 1 1 1 .41 to find a good F^-slice in each characteristic rectangle. 
Proposition 111.31 exports this good decay of the F|^-term to all constant w-slices 
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Figure 9: Closing the bootstrap. 



and hence the last outstanding bootstrap assumption ()132|) is finally retrieved with 
a better constant. 

What wc have shown is that A = A, so A is closed. This completes the proof of 
Proposition l7.ll The set A must therefore constitute the entire [0, oo] and hence the 
decay rates of Thcorcm l 1 . 1 1 are proven in the entire V, albeit in a different coordinate 
system than the one stated in the Theorem. The final subsection shows that the 
coordinate systems used in the bootstrap arc indeed close to the null-coordinate 
system defined in Theorem ll.il 



13.4 Convergence of Coordinate Systems 

What wc have already shown in section 18.3.21 is that the coordinates of a region 
^ ("i? (f,)) n {r* > rjj-} (a-priori defined in the coordinate system Cf,), are uniformly 
bounded in any coordinate system Cf for f > f. Iff] It is important to observe that 
the u-coordinatc in the region r* < is not uniformly close between the different 
coordinate systems. Indeed, in the coordinate system of Theorem 1 1 . 1 1 the horizon 
is located at w = oo, whereas in any coordinate system Cf it generically resides at 
a finite u value (eventually converging to m ^ oo for f — > oo). 

We finally establish the relation of the Cf to the coordinate system defined 
in Theorem 18.3.21 First recall that we have already shown that the geometrically 
defined point R of Theorem 11.11 (which features as an " origin" of the coordinate 
system) has coordinates uniformly close to {\/M, VM) in any coordinate system 
Cf , cf. section 18.31 In the second step we compare the scaling of the coordinates 
between the coordinate systems Cf and the one asserted by Theorem ll.il For this 
pick a point P on null-infinity. The value of 7 at this point in the coordinate system 
Cf can (for large enough f ) be estimated by integrating ([35|) from t = T (f) along 
a line of constant u: 

1<^(P)<\ + ^ (471) 

where is the area radius at the intersection of the Vr integral curve defining 
the coordinate system Cf and the null line u(P)@ In the limit f ^ 00 we have 
rjv ^ 00 and hence 7 (P) i. It follows that the scaling of the u-coordinate of Cf 
indeed converges to the one defined in Theorem 11.11 

The function k on the other hand satisfies |k — ^| < C (e) on 2? in both the 
coordinate systems Cf (cf . Proposition 14. 1|) and the one of Theorem 11.11 It is easy 
to show that with this bound holding on the null curve u — ^/M, the v coordinate 



**Note again that r = may change its location in the different coordinate systems but remains 
always close to the geometrically defined curve r = r/f of constant area radius. 
*^For this estimate only the smallness of 9 of Proposition 14. II is used. 
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of any two coordinate systems always satifies v ^ v for v > vM, which is aU what 
is needed to generahze decay statements in v to all coordinate systems. Namely 
integrating from the point W where the initial data intersect the null-line u = \fM 
(v ~ \fM there by previous remarks) to a point Q we have 

vq=vw+ dv < vw + C^ + C (e)) sup / r^ydv < -VM+A (tq ~ rw) 

Jvw u=\/II Jvw ^ 

and 

VQ=vw+ dv > vw+C^ - C {e)) inf / r.vdv>— — h- (rg - riy) . 

Jvw 'u=yA7 J- - M Jvw ^ ^ 

Hence w ~ w for any two coordinate systems. 

We have shown that the limit of the coordinate systems is a coordinate 
system in which the origin is slightly shifted compared to the one of Theorem 11.11 
and whose v scaling may be stretched or squeezed. It is now apparent that the 
decay rates stated also hold in the coordinate system of Theorem 11.11 



14 Final Comments and Open Questions 

Theorem 11.11 leaves room for generalizations. An obvious one is the treatment of 
the triaxial case, which at least conceptually is not expected to pose any difficulty. 
In fact the same vectorfields are expected to produce the required estimates for the 
fields B and C when contracted with an appropriate tensor T^^ - with the only 
additional catch coming from the coupling of B and C. A much more challenging 
problem is the derivation of better decay rates than the ones established here. As 
mentioned previously, in the context of compatible currents, the maximal decay rate 
is limited by the weights appearing in the if-vectorfield. It is an interesting question 
whether an additional vectorfield (or an entirely different idea) can extract stronger 
decay, which might be expected from the four-dimensional case [3]. An even more 
ambitious problem concerns the large data regime of the five-dimensional Bianchi 
IX model. The numerical studies of [T] suggest that a similar result to the one 
proven here should hold. In fact it may be possible to find an elaborate refinement 
of the ideas in [6j , which will allow an analysis of the large data regime within the 
symmetry class. 

Finally, there should exist various applications of the techniques to four-dimensional 
problems. As already mentioned in the introduction, the present paper may serve as 
a blueprint to obtain a small-data version of [6] for the self-gravitating scalar field. 
For genuinely novel results, the case of a conformally coupled scalar field could be 
investigated. 
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A Regularity and Green's identity 

It was remarked in section [3] that the coordinate systems Cf are . More precisely 
it was shown that they are piecewise with a discontinuity in spreading along 
the null-line v {B) and a discontinuity in along u (B). This discontinuity could 
be avoided by the introduction of a smooth interpolating function in the region 
around the cusp at the point B (cf. Figure [3]). However, as this would burden the 
notation even further, we will show here that the regularity is sufficient to carry out 
the calculations involving the vectorfields. 

Observing that the quantity P" defined in ([86|) is continuous and VqP" at least 
piecewise continuous (cf. ([92]) ). the basic identity (fTO|l is valid for the vectorfields 
X,Y,K in the coordinate systems Cf. 

For the vectorfields X and K we also make use of Green's identity (|105p in a 

region -"V\>^^. 




As depicted, the region may contain part of the null-line v {B) along which -j^ 
could be discontinuous and part of the null-line u (B) along which could be 
discontinuous. The functions D for which (|105p is applied are given by (|260p and 
(|373p . In both cases, D (u, v) is seen to be piecewise differentiable and such that OD 
is piecewise continuous. To derive the identity (|105p for these cases in our coordinate 

system, one should split the integration region ^[/^' t/j i^^to three pieces, along the 
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null lines u{B) and v{B), introducing additional boundary terms from the bold 
lines. Green's identity is then clearly valid in each subregion because all functions 
admit appropriate regularity there, i.e. in particular D is differentiable and DD is 
continuous in the interior. The integrand of the additional boundary term along 
the null-line v (B) however 

j [B^duD - DduB^] r^du (472) 

is continuous because the u derivative of D, which involves only the term (but 
not its w-analogue!), is continuous there. It is also bounded and the above integral 
will appear with a different sign for the two subregions. Analogously, the integrand 
of the other boundary term 

j [B^dyD ~ DdyB'^] r^dv (473) 

is continuous because the v derivative of D is continuous there. Hence adding up 
the three subregions the additional boundary terms cancel and the identity (jlOSp 
indeed holds as stated. 

B Different curves of constant r* 





very large and negative (close to the horizon). 




features as a source of smallness in the bootstrap 


<l 






negative, chosen in section [7.21 to make a certain bulk-term 




of the Y vectorfield positive in the region r* < 




functions a and (3 are supported in r* < —■^^/M only 


R* 


R* = -i^/M defined in Proposition |10.3| 


r* 

zero 


defined in section [HLZJ ~^VM < r*^^^ < -jqVM 





w 4Af (photon sphere for 5dim. Schwarzschild) 


R* 


squashing field on initial data is not supported for r* > R* 


R* 


defined in Lemma 112. Hequips a certain inteerand with a sien 




in a particular region 
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C Glossary 



a function depending on r*, used in the definition of the vectorficld Y 

P function depending on r*, used in the definition of the vectorficld Y 

B squashing field 

7 defined in ([33|) 

V defined in (gl]) 

5, 5 smallness parameters 

e, e smallness parameters 

C C = riB,, 

r] smallness parameter (cf. Corollarv l4.2l and Proposition [721) 

e e = r^B^^ 

d function used for the definition of the coordinate systems Cf, cf. (|43|) 

K defined in ([55]) 

A A = r^y 

M ^^ = T^ 

^ function depending on r* defined in (|277p 

m Hawking mass 

Mf final Bondi mass 

Ma Hawking mass at the point A, cf. section [3] 

r r (u, v) area radius 

p defined in ^ 

defined in dH) 

St defined in (fTTG]) 

(7 parameter, chosen in the section on the vectorfield X, cf. (|276p 

r affine parameter along = AM a, section [3] 

f affine parameter along = AAIf, section |3] 

, Lp2 defined in ([36)) 

X, X smooth interpolating functions, cf. (|137p and Proposition p0.3p 

ip smallness parameter, section [Y. 21 

il^ metric function 
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